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c^ I The interaction of an atom with an electromagnetic field is discussed in the 

presence of a time periodic external modulating force. It is explained that a 

control on atom by electromagnetic fields helps to design the quantum 

analog of classical optical systems. In these atom optical systems chaos 

may appear at the onset of external fields. The classical and quantum 

5^ I chaotic dynamics is discussed, in particular in an atom optics Fermi 

accelerator. It is found that the quantum dynamics exhibits dynamical 

localization and quantum recurrences. 



I Introduction 

Two hundred years ago, Gottingen physicist George Christoph Lichtenberg 
wrote "I think it is a sad situation in all our chemistry that we are unable 
to suspend the constituents of matter free" . Today, the possibilities to store 
atoms and to cool them to temperatures as low as micro kelvin and nano 
kelvin scale, have made atom optics a fascinating subject. It further provides 
a playground to study the newer effects of quantum coherence and quantum 
interference. 

In atom optics we take into account the internal and external degrees of 
freedom of an atom. The atom is considered as a de Broglie matter wave and 
these are optical fields which provide components, such as, mirrors, cavities 
and traps for the matter waves (Meystre 2001, Dowling and Gea-Banacloche 
1996). Thus, we find a beautiful manifestation of quantum duality. 

In atom optics systems another degree of freedom may be added by pro- 
viding an external periodic electromagnetic field. This arrangement makes 
it feasible to open the discussion on chaos. These periodically driven atom 
optics systems help to realize various dynamical systems which earlier were 
of theoretical interest only. 

The simplest periodically driven system, which has inspired the scien- 
tists over many decades to understand various natural phenomena, is Fermi 
accelerator. It is as simple as a ball bouncing elastically on a vibrating hori- 
zontal plane. The system also contributes enormously to the understanding 
of dynamical chaos (Lichtenberg 1980, 1983, 1992). 



I.l Atom optics: An overview 

In the last two decades, it has become feasible to perform experiments with 
cold atoms in the realm of atom optics. Such experiments have opened the 
way to find newer effects in the behavior of atoms at very low temperature. 
The recent developments in atom optics (Mlynek 1992, Adams 1994, Fil- 
let 1994, Arimondo 1996, Raizen 1999, Meystre 2001) make the subject a 
suitable framework to realize dynamical systems. 

Quantum duality, the work-horse of atom optics, is seen at work in ma- 
nipulating atoms. The atomic de Broglie waves are defiected, focused and 
trapped. However, these are optical fields which provide tools to manipulate 
the matter waves. As a consequence, in complete analogy to classical optics, 
we have atom optical elements for atoms, such as, mirrors, beam splitters, 
interferometers, lenses, and waveguides (Sigel 1993, Adams 1994, Dowling 
1996, Theuer 1999, Meystre 2001). 

As a manifestation of the wave-particle duality, we note that a standing 
light wave provides optical crystal (Sleator 1992a, Sleator 1992b). Thus, we 
may find Raman-Nath scattering and Bragg scattering of matter waves from 
an optical crystal (Saif 2001a, Khalique 2003). In addition, an exponen- 
tially decaying electromagnetic field acts as an atomic mirror (Balykin 1988, 
Kasevich 1990, Wallis 1992). 

Atom interferometry is performed as an atomic de Broglie wave scatters 
through two standing waves acting as optical crystal, and aligned parallel 



to each other. The matter wave sphts into coherent beams which later re- 
combine and create an atom interferometer (Rasel 1995). The atomic phase 
interferometry is performed as an atomic de Broghe wave reflects back from 
two different positions of an atomic mirror, recombines, and thus interferes 
(Henkel 1994, Steane 1995, Szriftgiser 1996). 

An atom optical mirror for atoms (Kazantsev 1974, Kazantsev 1990) is 
achieved by the total internal reflection of laser light in a dielectric slab 
(Balykin 1987, Balykin 1988). This creates, an exponentially decaying elec- 
tromagnetic field appears outside of the dielectric surface. The decaying field 
provides an exponentially increasing repulsive force to a blue detuned atom, 
which moves towards the dielectric. The atom exhausts its kinetic energy 
against the optical field and reflects back. 

For an atom, which moves under gravity towards an atomic mirror, the 
gravitational field and the atomic mirror together act like a cavity — named 
as an atomic trampoline (Kasevich 1990) or a gravitational cavity (Wallis 
1992). The atom undergoes a bounded motion in the system. 

It is suggested by H. Wallis that a small change in the curvature of the 
atomic mirror helps to make a simple surface trap for the bouncing atom 
(Wallis 1992). An atomic mirror, comprising a blue detuned and a red de- 
tuned optical field with different decay constants, leads to the atomic trap- 
ping as well (Ovchinnikov 1991, DesbioUes 1996). 

The experimental observation of the trapping of atoms over an atomic 
mirror in the presence of gravitational field was made by ENS group in Paris 



Figure 1: Observation of trapping of atoms in a gravitational cavity: (a) 
Schematic diagram of the experimental set up. (b) Number of atoms detected 
by the probe beam after their initial release as a function of time are shown 
as dots. The solid curve is a result of corresponding Monte-Carlo simulations. 
(Aminoff 1993). 

(Aminoff 1993). In the experiment cold cesium atoms were dropped from a 
magneto-optic trap on an atomic mirror, developed on a concave spherical 
substrate, from a height of 3mm. The bouncing atoms were observed more 
than eight times, as shown in figure 1. 

A fascinating achievement of the gravitational cavity is the development 
of recurrence tracking microscope (RTM) to study surface structures with 
nano- and sub nano-meter resolutions (Saif 2001). The microscope is based 
on the phenomena of quantum revivals. 

In RTM, atoms observe successive reflections from the atomic mirror. The 
mirror is joined to a cantilever which has its tip on the surface under investi- 
gation. As the cantilever varies its position following the surface structures, 
the atomic mirror changes its position in the upward or downward direc- 
tion. The time of a quantum revival depends upon the initial height of the 
atoms above the mirror which, thus, varies as the cantilever position changes. 
Hence, the change in the time of revival reveals the surface structures under 
investigation. 

The gravitational cavity has been proposed (Ovchinnikov 1995, Soding 



1995, Laryushin 1997, Ovchinnikov 1997) to cool atoms down to the micro 
Kelvin temperature regime as well. Further cooling of atoms has made it 
possible to obtain Bose-Einstein condensation (Davis 1995, Anderson 1995, 
Bradley 1995), a few micrometers above the evanescent wave atomic mirror 
(Hammes 2003, Rychtarik 2004). 

A modulated gravitational cavity constitutes atom optics Fermi acceler- 
ator for cold atoms(Chen 1997, Saif 1998). The system serves as a suitable 
framework to analyze the classical dynamics and the quantum dynamics in 
laboratory experiments. A bouncing atom displays a rich dynamical behavior 
in the Fermi accelerator (Saif 1999). 

1.2 The Fermi accelerator 

In 1949, Enrico Fermi proposed a mechanism for the mysterious origin of 
acceleration of the cosmic rays (Fermi 1949). He suggested that it is the 
process of collisions with intra-galactic giant moving magnetic fields that 
accelerates cosmic rays. 

The accelerators based on the original idea of Enrico Fermi display rich 
dynamical behavior both in the classical and the quantum evolution. In 1961, 
Ulam studied the classical dynamics of a particle bouncing on a surface which 
oscillates with a certain periodicity. The dynamics of the bouncing particle 
is bounded by a fixed surface placed parallel to the oscillating surface (Ulam 
1961). 



In Fermi-Ulam accelerator model, the presence of classical chaos (Lieber- 
man 1972, Lichtenberg 1980, 1983, 1992) and quantum chaos (Karner 1989, 
Seba 1990) has been proved. A comprehensive work has been devoted to 
study the classical and quantum characteristics of the system (Lin 1988, 
Makowski 1991, Reichl 1992, Dembihski 1995). However, a particle bounc- 
ing in this system has a limitation, that, it does not accelerate forever. 

Thirty years after the first suggestion of Fermi, Pustyl'nikov provided 
detailed study of another accelerator model. He replaced the fixed horizontal 
surface of Fermi-Ulam model by a gravitational field. Thus, Pustyl'nikov 
considered the dynamics of a particle on a periodically oscillating surface in 
the presence of a gravitational field (Pustyl'nikov 1977). 

In his work, he proved that a particle bouncing in the accelerator system 
finds modes, where it ever gets unbounded acceleration. This feature makes 
the Fermi-Pustyl'nikov model richer in dynamical beauties. The schematic 
diagram of Fermi-Ulam and Fermi-Pustyl'nikov model is shown in Fig 2. 

In case of the Fermi-Ulam model, the absence of periodic oscillations of 
reflecting surface makes it equivalent to a particle bouncing between two 
fixed surfaces. However, in case of the Fermi-Pustyl'nikov model, it makes 
the system equivalent to a ball bouncing on a fixed surface under the in- 
fiuence of gravity. These simple systems have thoroughly been investigated 
in classical and quantum domains (Langhoff 1971, Gibbs 1975, Desko 1983, 
Goodins 1991, Whineray 1992, Seifert 1994, Bordo 1997, Andrews 1998, Gea- 
Banacloche 1999). 
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Figure 2: a) Schematic diagram of the Fermi-Ulam accelerator model: A 
particle moves towards a periodically oscillating horizontal surface, experi- 
ences an elastic collision, and bounces off in the vertical direction. Later, 
it bounces back due to another fixed surface parallel to the previous one. 
b) Schematic diagram of the Fermi- Pustyl'nikov model: A particle observes 
a bounded dynamics on a periodically oscillating horizontal surface in the 
presence of a constant gravitational field, g. The function, f(t), describes the 
periodic oscillation of the horizontal surfaces. 

In the presence of an external periodic oscillation of the refiecting surface, 
the Fermi-Ulam model and the Fermi- Pustyl'nikov model display the mini- 
mum requirement for a system to be chaotic (Lichtenberg 1983, 1992). For 
the reason, these systems set the stage to understand the basic characteristics 
of the classical and quantum chaos (Jose 1986, Seba 1990, Badrinarayanan 
1995, Mehta 1990, Reichl 1986). Here, we focus our attention mainly on the 
classical and quantum dynamics in the Fermi-Pustyl'nikov model. For the 
reason, we name it as Fermi accelerator model in the rest of the report. 

1.3 Classical and quantum chaos in atom optics 

Quantum chaos, as the study of the quantum characteristics of the classically 
chaotic systems, got immense attention after the work of Bayfield and Koch 
on microwave ionization of hydrogen (Bayfield and Koch, 1974). In the 
system the suppression of ionization due to the microwave field was attributed 
to dynamical localization (Casati et al.,1987;Koch and Van Leeuwen 1995). 



Later, the phenomenon was observed experimentally (Galvez 1988, Bliimel 
1989, Bayfield 1989, Arndt 1991, Benvenuto 1997, Segev 1997). 

Historically, the pioneering work of Giulio Casati and co-workers on delta 
kicked rotor unearthed the remarkable property of dynamical localization of 
quantum chaos (Casati 1979). They predicted that a quantum particle ex- 
hibits diffusion following classical evolution till quantum break time. Beyond 
this time the diffusion stops due to quantum interference effects. In 1982, 
Fishman, Grempel and Prange proved mathematically that the phenomenon 
of the dynamical localization in kicked rotor is the same as Anderson local- 
ization of solid state physics (Fishman 1982). 

The study of quantum chaos in atom optics began with a proposal by 
Graham, Schlautmann, and Zoller (Graham 1992). They investigated the 
quantum characteristics of an atom which passes through a phase modulated 
standing light wave. During its passage the atom experiences a momentum 
transfer by the light field. The classical evolution in the system exhibits 
chaos and the atom displays diffusion. However, in the quantum domain, 
the momentum distribution of the atom at the exit is exponentially localized 
or dynamically localized, as shown in figure 3. 

Experimental study of quantum chaos in atom optics is largely based 
upon the work of Mark Raizen (Raizen 1999). In a series of experiments he 
investigated the theoretical predictions regarding the atomic dynamics in a 
modulated standing wave field and regarding delta kicked rotor model. The 



Figure 3: Dynamical localization of an atom in a phase modulated standing 
wave field: The time averaged probability distribution in momentum space 
displays the exponentially localized nature of momentum states. The ver- 
tical dashed lines give the border of the classically chaotic domain. The 
exponential behavior on semi-logarithmic plot is shown also by the dashed 
lines. (Graham et al. 1992). 

work also led to the invention of newer methods of atomic cooling as well 
(Ammann 1997). 

In the framework of atom optics, periodically driven systems have been 
explored to study the characteristics of the classical and quantum chaos. 
These systems include an atom in a modulated electro-magnetic standing 
wave field (Graham 1992, Raizen 1999), an ion in a Paul trap in the presence 
of an electromagnetic field (Ghafar 1997), an atom under the infiuence of 
strong electromagnetic pulses (Raizen 1999), and an atom in a Fermi accel- 
erator (Saif 1998, 2000, 2000a, 2002). 

The atom optics Fermi accelerator is advantageous in many ways: It is 
analogous to the problem of a hydrogen atom in a microwave field. In the 
absence of external modulation, both the systems possess weakly binding 
potentials for which level spacing reduces with increase in energy. However, 
Fermi accelerator is more promising due to the absence of inherent continuum 
in the unmodulated case as it is found in the hydrogen atom. 

For a small modulation strength and in the presence of a low frequency 
of the modulation, an atom exhibits bounded and integrable motion in the 
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classical and quantum domain (Saif 1998). However, for larger values of the 
strength and/or the higher frequency of the modulation, there occurs clas- 
sical diffusion. In the corresponding quantum domain, an atom displays no 
diffusion in the Fermi accelerator, and eventually displays exponential local- 
ization both in the position space and the momentum space. The situation 
prevails till a critical value of the modulation strength which is based purely 
on quantum laws. 

The quantum delocalization (Chirikov 1986) of the matter waves in the 
Fermi accelerator occurs at higher values of strength and frequency of the 
modulation, above the critical value. The transition from dynamical local- 
ization to quantum delocalization takes place as the spectrum of the Floquet 
operator displays transition from pure point to quasi-continuum spectrum 
(Brenner 1996, Oliveira 1994, Benvenuto 1991). 

In nature, interference phenomena lead to revivals (Averbukh 1989a, 
Averbukh 1989b, Alber 1990, Fleischhauer 1993, Chen 1995, Leichtle 1996a, 
Leichtle 1996b, Robinett 2000). The occurrence of the revival phenomena in 
time dependent systems has been proved to be their generic property (Saif 
2005), and regarded as a test of deterministic chaos in quantum domain 
(Bliimel 1997). The atomic evolution in Fermi accelerator displays revival 
phenomena as a function of the modulation strength and initial energy of 
the atom (Saif 2000, 2000c, 2000d, 2000e, 2002). 
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1.4 Layout 

This report is organized as follows: In section II, we review the interac- 
tion of an atom with an optical field. In section III, we briefiy summarize 
the essential ideas of the quantum chaos dealt within the report. The experi- 
mental progress in developing atom optics elements for the atomic de Broglie 
waves is discussed in section IV. These elements are the crucial ingredients 
of complex atom optics systems, as explained in section V. We discuss the 
spatio-temporal characteristics of the classical chaotic evolution in section VI. 
The Fermi accelerator model is considered as the focus of our study. In the 
corresponding quantum system the phenomenon of dynamical localization is 
discussed in section VII. The study of the recurrence phenomena in general 
periodically driven systems is presented in section IX, and their study in the 
Fermi accelerator is made in section X. In section XI, we make a discussion 
on chaotic dynamics in periodically driven atom optics systems. We conclude 
the report, in section XII, by a brief discussion of decoherence in quantum 
chaotic systems. 

II Interaction of an atom with an optical field 

In this section, we provide a review of the steps leading to the effective 
Hamiltonian which governs atom-field interaction. We are interested in the 
interaction of a strongly detuned atom with a classical light field. For the 
reason, we provide quantum mechanical treatment to the center of mass 
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motion however develop semi-classical treatment for the interaction of the 
internal degrees of freedom and optical field. 

For a general discussion of the atom-field interaction, we refer to the lit- 
erature (Cohen- Tannoudji 1977, Kazantsev 1990, Meystre 2001, Shore 1990, 
Sargent 1993, Scully 1997, Schleich 2001). 

II. 1 Interacting atom 

We consider a two-level atom with the ground and the excited states as \g) 
and |e), respectively. These states are eigen states of the atomic Hamiltonian, 
Hq, and the corresponding eigen energies are huj^g-^ and huj(^e)- Thus, in the 
presence of completeness relation, \g){g\ + |e)(e| = 1, we may define the 
atomic Hamiltonian, as 

Ho = hu(e)\e){e\+hLj(^g)\g){g\. (1) 

In the presence of an external electromagnetic field, a dipole is formed 
between the electron, at position re, and the nucleus of the atom at, r„. 
In relative coordinates the dipole moment becomes er^ = e(re — r„). Conse- 
quently, the interaction of the atom with the electromagnetic field is governed 
by the interaction Hamiltonian, 

H,nt = -ero-E{r + 6r,t). (2) 

Here, r is the center-of-mass position vector and Sr is either given by Sr = 
—iTLeV / {iTLe + iTin) or 5y = iTinTc / {iTie + mn) ■ The symbols nie and m„ define 
mass of the electron and that of nucleus, respectively. 
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II. 1.1 Dipole Approximation 

Keeping in view a comparison between tlie wavelength of the electromagnetic 
field and the atomic size, we consider that the field does not change signifi- 
cantly over the dimension of the atom and take E{r + 6r, t) = E{r, t). Thus, 
the interaction Hamiltonian becomes, 

Hint^ -evo- E{v,t). (3) 

In the analysis we treat the atom quantum mechanically, hence , we 
express the dipole in operator description, as p. With the help of the com- 
pleteness relation, we may define the dipole operator, as 

^=1.^.1 = e{\g){g\ + |e)(e|)r„(|(?)((?| + |e)(e|) = p\e){g\ + p'\g){e\. (4) 

Here, we have introduced the matrix element, p = e(e|ro|5'), of the electric 
dipole moment. Moreover, we have used the fact that the diagonal elements, 
{g\ro\g) and (e|ro|e), vanish as the energy eigen states have well-defined 
parity. 

We define a^ = \e){g\ and a = \g){e\, as the atomic raising and lowering 
operators, respectively. In addition, we may consider the phases of the states 
|e) and \g) such that the matrix element p is real (Kazantsev 1990). This 
leads us to represent the dipole operator, as 

p = p{a + a^). (5) 
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Therefore, we arrive at the interaction Hamiltonian 



i^mt = -((xt + <T)p-E(r,t). (6) 

Thus, the Hamihonian which descries the interaction of the atom with 
the optical field in the presence of its center-of-mass motion, reads as 

H=^ + %e)|e)(e| + hLU^g)\g){g\ - {a^ + a)p- E(f,t), (7) 

where, p describes the center-of-mass momentum. Here, the first term on the 
right hand side corresponds to the kinetic energy. It becomes crucial when 
there is a significant atomic motion during the interaction of the atom with 
the optical field. 

1 1. 2 Effective potential 

The evolution of the atom in the electromagnetic field becomes time depen- 
dent in case the field changes in space and time, such that 

E(r, t) = Eq u(r) cos Uft. (8) 

Here, Eq and uof are the amplitude and the frequency of the field, respectively. 
The quantity u(r) = eo'u(r) denotes the mode function, where Bq describes 
the polarization vector of the field. 

II. 2.1 Rotating wave approximation 

The time dependent Schrodinger equation, 

*M . HW, (9) 
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controls the evolution of the atom in the time dependent electromagnetic 
field. Here, H, denotes the total Hamiltonian, given in equation (7), in the 
presence of field defined in equation (8). 

In order to solve the time dependent Schrodinger equation, we express 
the wave function, \%l)), as an ansatz, i.e. 



g-u 



^^^%{r,t)\g) + e-^H«)+-/)Ve(r,t)|e). (10) 



Here, ipg = ipg{r,t) and ipe = i^ei^'yt), express the probability amplitudes in 
the ground state and excited state, respectively. 

We substitute equation (10), in the time dependent Schrodinger equation, 
given in equation (9). After a little operator algebra we, finally get the 
coupled equations for the probability amplitudes, ipg and ipe, expressed as 

'^^ = ^^9 - l^^Ru{r){l + e-2^/*)^e, (11) 

^h^ = l^^e - Mt^e - lnnnu{r){l + e~'^f')ijg. (12) 

Here, ^2^ = p ■ eQEQ/2h, denotes the Rabi frequency, and S = Uf — Ueg is 
the measure of the tuning of the external field frequency, Uf, away from the 
atomic transition frequency, expressed as Ueg = u;(e) — uj(^g). 

In the rotating wave approximation, we eliminate the rapidly oscillating 
terms in the coupled equations (11) and (12) by averaging them, over a period, 
T = 2'K/ujf. We assume that the probability amplitudes, ipe and ipg, do not 
change appreciably over the time scale and approximate, 

l + e-2^^*~l. (13) 
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Thus, we eliminate the exphcit time dependence in the equations (11) and 
(12), and get 

^h^ = ^^g - hnnu{r)^e, (14) 

Oil) p 
z/i— ^ = — ^e - fiSipe - hQRu{r)iljg. (15) 

II. 2. 2 Adiabatic approximation 

In the hmit of a large detuning between the field frequency, Uf, and the 
atomic transition frequency uJeg, we consider that the excited state population 
changes adiabatically. Hence, we take ipe{t) ~ V^e(O) = constant. 

As a consequence, in equation (15) the deviation of the probability am- 
plitude ipe with respect to time and position becomes vanishingly small. This 
provides an approximate value for the probability amplitude ipeit), as 

Mt) ^ -^u{r)^g{t). (16) 

Thus, the time dependent Schrodinger equation for the ground state prob- 
ability amplitude, ipg, becomes 

For simplicity, here and later in the report, we drop subscript g and take 

The Schrodinger equation, given in equation (17), effectively governs the 
dynamics of a super-cold atom, in the presence of a time dependent optical 
field. The atom almost stays in its ground state under the condition of a 
large atom-field detuning. 
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II.2.3 The Effective Hamiltonian 

Equation (17) leads us to an effective Hamiltonian, -ffg//; such as 

Here, the first term describes the center-of-mass kinetic energy, whereas, 
the second term indicates an effective potential as seen by the atom. The 
spatial variation in the potential enters through the mode function of the 
electromagnetic field, u. Fascinatingly, we note that by a proper choice of 
the mode function almost any potential for the atom can be made. 

The second term in equation (18), leads to an effective force F = —2{hVt^ / 5)u\/ 
u, where \/u describes gradient of u. Hence, the interaction of the atom with 
the electromagnetic field exerts a position dependent gradient force on the 
atom, which is directly proportional to the square of the Rabi frequency, ^2^, 
and inversely proportional to the detuning, S. 

II. 3 Scattering atom 

In the preceding subsections we have shown that an atom in a spatially vary- 
ing electromagnetic field experiences a position dependent effective potential. 
As a consequence the atom with a dipole moment gets deflected as it moves 
through an optical field (Moskowitz 1983) and experiences a position depen- 
dent gradient force. The force is the largest where the gradient is the largest. 
For example, an atom in a standing light field experiences a maximum force 
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at the nodes and at antinodes it observes a vanishing dipole force as the 
gradient vanishes. 

We enter into a new regime, when the kinetic energy of the atom paral- 
lel to the optical field is comparable to the recoil energy. We consider the 
propagation of the atom in a potential formed by the interaction between 
the dipole and the field, moreover, we treat the center-of-mass motion along 
the cavity quantum mechanically. 

Ill Quantum characteristics of chaos 

Over last twenty five years the field of atom optics has matured and become a 
rapidly growing branch of quantum electronics. The study of atomic dynam- 
ics in a phase modulated standing wave by Graham, Schlautman and Zoller 
(Graham 1992), made atom optics a testing ground for quantum chaos as 
well. Their work got experimental verification as the dynamical localization 
of cold atoms was observed later in the system in momentum space (Moore 
1994, Bardroff 1995, Amman 1997). 

III.l Dynamical localization 

In 1958, P. W. Anderson showed the absence of diffusion in certain random 
lattices (Anderson 1958). His distinguished work was recognized by a Nobel 
prize, in 1977, and gave birth to the phenomenon of localization in solid state 
physics — appropriately named as Anderson localization (Anderson 1959). 



An electron, on a one dimensional crystal lattice displays localization if 
the equally spaced lattice sites are taken as random. The randomness may 
arise due to the presence of impurities in the crystal. Thus, at each ith site of 
the lattice, there acts a random potential, Tj. The probability amplitude of 
the hopping electron, Ui, is therefore expressed by the Schrodinger equation 

TiUi + Y^ WrUi+r = EUi, (19) 

where, Wr is the hopping amplitude from ith site to its rth neighbor. If the 
potential, Ti, is periodic along the lattice, the solution to the equation (19) 
is Bloch function, and the energy eigen values form a sequence of continuum 
bands. However, in case Ti are uncorrelated from site to site, and distributed 
following a distribution function, eigen values of the equation (19) are expo- 
nentially localized. Thus, in this situation the hopping electron finds itself 
localized. 

Twenty years later, Giulio Casati and coworkers suggested the presence of 
a similar phenomenon occurring in the kicked rotor model (Casati 1979). In 
their seminal work they predicted that a quantum particle, subject to peri- 
odic kicks of varying strengths, exhibits the suppression of classical diffusion. 
The phenomenon was named as dynamical localization. Later, on mathemat- 
ical grounds, Fishman, Grempel and Prange developed equivalence between 
dynamical localization and the Anderson localization (Fishman 1982). 

In its classical evolution an ensemble of particles displays diffusion as it 
evolves in the kicked rotor in time. In contrast, the corresponding quantum 
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system follows classical evolution only up to a certain time. Later, it dis- 
plays dynamical localization as the suppression of the diffusion or its strong 
reduction (Brivio 1988, Izrailev 1990). 

The maximum time for which the quantum dynamics mimics the cor- 
responding classical dynamics is the quantum break time of the dynamical 
system (Haake 1992, Bliimel 1997). Before the quantum break time, the 
quantum system follows classical evolution and fills the stable regions of the 
phase space (Reichl 1986). 

The phenomenon of dynamical localization is a fragile effect of coherence 
and interference. The quantum suppression of classical diffusion does not oc- 
cur under the condition of resonance or in the presence of some translational 
invariance (Lima 1991, Guarneri 1993). The suppression of diffusion in the 
absence of these conditions is dominantly because of restrictions of quantum 
dynamics by quantum cantori (Casati 1999a, Casati 1999b). 

The occurrence of dynamical localization is attributed to the change in 
statistical properties of the spectrum (Berry 1981, Zaslavsky 1981, Haake 
1992, Altland 1996). The basic idea involved is that, integrability cor- 
responds to Poisson statistics (Brody 1981, Bohigas 1984), however, non- 
integrability corresponds to Gaussian orthogonal ensemble (GOE) statistics 
as a consequence of the Wigner's level repulsion. 

In the Fermi-Ulam accelerator model, for example, the quasi energy spec- 
trum of the Floquet operator displays such a transition. It changes from Pois- 
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son statistics to GOE statistics as the effective Planck's constant changes in 
value (Jose 1986). 

The study of the spectrum leads to another interesting understanding of 
localization phenomenon. Under the influence of external periodic force, as 
the system exhibits dynamical localization the spectrum changes to a pure 
point spectrum (Prange 1991, Jose 1991, Dana 1995). There occurs a phase 
transition to a quasi continuum spectrum leading to quantum delocalization 
(Benvenuto 1991, Oliveira 1994), for example, in the Fermi accelerator. 

III. 2 Dynamical recurrences 

In one dimensional bounded quantum systems the phenomena of quantum 
recurrences exist (Robinett 2000). The roots of the phenomena are in the 
quantization or discreteness of the energy spectrum. Therefore, their occur- 
rence provides a profound manifestation of quantum interference. 

In the presence of an external time dependent periodic modulation, these 
systems still exhibit the quantum recurrence phenomena (Saif 2002, Saif 
2005). Interestingly, their occurrence is regarded as a manifestation of de- 
terministic quantum chaos (Bliimel 1997). 

Hogg and Huberman provided the first numerical study of quantum recur- 
rences in chaotic systems by analyzing the kicked rotor model (Hogg 1982). 
Later, the recurrence phenomena was further investigated in various clas- 
sically chaotic systems, such as, the kicked top model (Haake 1992), in the 
dynamics of a trapped ion interacting with a sequence of standing wave pulses 
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(Breslin 1997), in the stadium billiards (Tomsovic 1997) and in multi-atomic 
molecules (Grebenshchikov 1997). 

The phenomena of quantum recurrences are established generic to peri- 
odically driven quantum systems which may exhibit chaos in their classical 
domain (Saif 2002, Saif 2005). The times of quantum recurrences are cal- 
culated by secular perturbation theory, which helps to understand them as 
a function of the modulation strength, initial energy of the atom, and other 
parameters of the system (Saif 2000b, Saif 2000c). 

In a periodically driven system, the occurrence of recurrences depends 
upon the initial conditions in the phase space. This helps to analyze quan- 
tum nonlinear resonances and the quantum stochastic regions by studying 
the recurrence structures (Saif 2000c). Moreover, it is suggested that the 
quantum recurrences serve as a very useful probe to analyze the spectrum of 
the dynamical systems (Saif 2005). 

In the periodically driven one-dimensional systems, the spectra of Floquet 
operators serve as quasi-energy spectra of the time dependent systems. A 
lot of mathematical work has been devoted to the study of Floquet spectra 
of periodically driven systems (Breuer 1989, 1991a, 1991b). 

It has been established that, in the presence of external modulation in 
the tightly binding potentials where the level spacing between adjacent levels 
increases with the increase in energy, the quasi energy spectrum remains a 
point spectrum regardless of the strength of the external modulation (Haw- 
land 1979, 1987, 1989a, 1989b, Joye 1994). However, in the presence of 
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external modulation in weakly binding potentials, where the level spacing 
decreases with increase in energy it changes from a point spectrum to a con- 
tinuum spectrum above a certain critical modulation strength (Delyon 1985, 
Benvenuto 1991, Oliveira 1994, Brenner and Fishman 1996). In tightly bind- 
ing potentials, the survival of point spectrum reveals quantum recurrences, 
whereas, the disappearance of point spectrum in weakly binding potentials 
leads to the disappearance of quantum recurrences above a critical modula- 
tion strength (Saif 2000c, Iqbal 2005). 

III. 3 Poincare' recurrences 

According to the Poincare' theorem a trajectory always return to a region 
around its origin, but the statistical distribution depends on the dynamics. 
For a strongly chaotic motion, the probability to return or the probability 
to survive decays exponentially with time (Lichtenberg 1992). However, in 
a system with integrable and chaotic components the survival probability 
decays algebraically (Chirikov 1999). This change in decay rate is subject to 
the slowing down of the diffusion due to chaos border determined by some 
invariant tori (Karney 1983, Chirikov 1984, Meiss 1985, Meiss 1986, Chirikov 
1988, Ruffo 1996). 

The quantum effects modify the decay rate of the survival probability 
(Tanabe 2002). It is reported that in a system with phase space a mixture 
of integrable and chaotic components, the algebraic decay P{t) ~ 1/t^ has 
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the exponent p = 1. This behavior is suggested to be due to tunnehng and 
locahzation effects (Casati 1999c). 

IV Mirrors and cavities for atomic de Broglie 
waves 

As discussed in section II, by properly tailoring the spatial distribution of 
an electromagnetic field we can create almost any potential we desire for 
the atoms. Moreover, we can make the potential repulsive or attractive by 
making a suitable choice of the atom-field detuning. An atom, therefore, 
experiences a repulsive force as it interacts with a blue detuned optical field, 
for which the field frequency is larger than the transition frequency. However, 
there is an attractive force on the atom if it finds a red detuned optical field, 
which has a frequency smaller than the atomic transition frequency. Thus, 
in principle, we can construct any atom optics component and apparatus for 
the matter waves, analogous to the classical optics. 

IV. 1 Atomic mirror 

A mirror for the atomic de Broglie waves is a crucial ingredient of the atomic 
cavities. An atomic mirror is obtained by an exponentially decaying optical 
field or an evanescent wave field. Such an optical field exerts an exponen- 
tially increasing repulsive force on an approaching atom, detuned to the blue 
(Bordo 1997). 
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How to generate an evanescent wave field is indeed an interesting ques- 
tion. In order to answer this question, we consider an electromagnetic field 
E{r,t) = £'(r)e~*'^^*, which travels in a dielectric medium with a dielectric 
constant, n and undergoes total internal refiection. The electromagnetic field 
inside the dielectric medium reads 

E{r,t) = Eoe'^-'""'f'er. (20) 

where, e,. is the polarization vector and k = kk is the propagation vector. 

The electromagnetic field, E{r,t), is incident on an interface between 
the dielectric medium with the dielectric constant, n, and another dielectric 
medium with a smaller dielectric constant, rii. The angle of incidence of the 
field is 6i with the normal to the interface. 

Since the index of refraction rii is smaller than n, the angle Oj. at which 
the field refracts in the second medium, is larger than 6i. As we increase 
the angle of incidence 6i, we may reach a critical angle, 6i = 6^, for which 
6r = 7i/2. According to Snell's law, we define this critical angle of incidence 
as 

e, ^ sin- (^) . (21) 

Hence, for an electromagnetic wave with an angle of incidence larger than 
the critical angle, that is, 6i > 6c, we find the inequality, sin^,. > 1 (Mandel 
1986, Mandel and Wolf 1995). As a result, we deduce that 9r is imaginary, 
and define 

cos Or = i^ 



( sin Oi \ 



\ Vsin6'c 
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1. (22) 



Therefore, the field in the medium of smaller refractive index, rii, reads 

E{r, t) = e^Eoe'^''-''"''^'-+'^''^"''^'-e-''^f' = erEoe"''^e*(^^-'^^*) (23) 



where, k = /ciW (sin^j/sin^c)^ — 1 and P = kisinOi/ sinOc (Jackson 1965). 
Here, ki defines the wave number in the medium with the refractive index 

III. This demonstrates that in case of total internal reflection the field along 
the normal of the interface decays in the positive z direction, in the medium 
with the smaller refractive index. 

In 1987, V. Balykin and his coworkers achieved the first experimental 
realization of an atomic mirror (Balykin 1987). They used an atomic beam 
of sodium atoms incident on a parallel face plate of fused quartz and observed 
the specular refiection. 

They showed that at small glancing angles, the atomic mirror has a re- 
flection coefficient equal to unity. As the incident angle increases a larger 
number of atoms reaches the surface and undergoes diffusion. As a result 
the reflection coefficient decreases. The refiection of atoms bouncing perpen- 
dicular to the mirror is investigated in reference (Aminoff 1993) and from a 
rough atomic mirror studied in reference (Henkel 1997). 

IV. 1.1 Magnetic mirror 

We can also construct atomic mirror by using magnetic fields instead of 
optical fields. At first, magnetic mirror was realized to study the reflection 
of neutrons (Vladimirskii 1961). In atom optics, the use of a magnetic mirror 
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was suggested in reference (Opat 1992), and later it was used to study the 
reflection of incident rubidium atoms perpendicular to the reflecting surface 
(Roach 1995, Hughes 1997a, Hughes 1997b, Saba 1999). Recently, it has 
become possible to modulate the magnetic mirror by adding a time dependent 
external field. We may also make controllable corrugations which can be 
varied in a time shorter than the time taken by atoms to interact with the 
mirror (Rosenbusch 2000a, Rosenbusch 2000b). 

Possible mirror for atoms is achieved by means of surface plasmons as well 
(Esslinger 1993, Feron 1993, Christ 1994). Surface plasmons are electromag- 
netic charge density waves propagating along a metallic surface. Traveling 
light waves can excite surface plasmons. The technique provides a tremen- 
dous enhancement in the evanescent wave decay length (Esslinger 1993). 

IV. 2 Atomic cavities 

Based on the atomic mirror various kinds of atomic cavities have been sug- 
gested. A system of two atomic mirrors placed at a distance with their ex- 
ponentially decaying fields in front of each other form a cavity or resonator 
for the de Broglie waves. The atomic cavity is regarded as an analog of the 
Fabry Perot cavity for radiation fields (Svelto 1998). 

By using more than two mirrors, other possible cavities can be developed 
as well. For example, we may create a ring cavity for the matter waves by 
combining three atomic mirrors (Balykin 1989). 
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An atomic gravitational cavity is a special arrangement. Here, atoms 
move under gravity towards an atomic mirror, made up of an evanescent 
wave (Matsudo 1997). The mirror is placed perpendicular to the gravita- 
tional field. Therefore, the atoms observe a normal incidence with the mirror 
and bounce back. Later, they exhaust their kinetic energy moving against 
the gravitational field (Kasevich 1990, Wallis 1992) and return. As a con- 
sequence, the atoms undergo a bounded motion in this atomic trampoline 
or atomic gravitational cavity. Hence, the evanescent wave mirror together 
with the gravitational field constitutes a cavity for atoms. 

IV. 3 Gravitational cavity 

The dynamics of atoms in the atomic trampoline or atomic gravitational 
cavity attracted immense attention after the early experiments by the group 
of S. Chu at Stanford, California (Kasevich 1990). They used a cloud of 
sodium atoms stored in a magneto-optic trap and cooled down to 25;uK. 

As the trap switched off the atoms approached the mirror under grav- 
ity and display a normal incidence. In their experiments two bounces of 
the atoms were reported. The major noise sources were fiuctuations in the 
laser intensities and the number of initially stored atoms. Later experiments 
reported up to thousand bounces (Ovchinnikov 1997). 

Another kind of gravitational cavity can be realized by replacing the 
optical evanescent wave field by liquid helium, forming the atomic mirror for 
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hydrogen atoms. In this setup hydrogen atoms are cooled below 0.5K° and 
a specular reflection of 80% has been observed (Berkhout 1989). 

IV. 3.1 A bouncing atom 

An atom dropped from a certain initial height above an atomic mirror expe- 
riences a linear gravitational potential, 

Vgr = mgz, (24) 

as it approaches the mirror. Here, m denotes mass of the atom, g expresses 
the constant gravitational acceleration and z describes the atomic position 
above the mirror. Therefore, by taking the optical field given in Eq. 23 and 
using Dipole approximation and rotating wave approximation the effective 
Hamiltonian, given in equation (18), becomes 

Heff = -^— + mgz + ^fiez/e"^"^. (25) 

^ fit 

The effective Hamiltonian governs the center-of-mass motion of the atom in 
the presence of the gravity above evanescent wave field. Here, ^e// = ^1?/^) 
describes the effective Rabi frequency, and /t^^ describes the decay length of 
the atomic mirror. 

The optical potential is dominant for smaller values of z and decays ex- 
ponentially as z becomes larger. Thus, for larger positive values of z as the 
optical potential approaches zero the gravitational potential takes over, as 
shown in figure 4 (left). 
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Figure 4: (left) As we switch off the Magneto-optic trap (MOT) at time, 
t = 0, the atomic wave packet starts its motion from an initial height. It 
moves under the influence of the linear gravitational potential Vgj. (dashed 
line) towards the evanescent wave atomic mirror. Close to the surface of the 
mirror the effect of the evanescent light field is dominant and the atom expe- 
riences an exponential repulsive optical potential V^pt (dashed line). Both the 
potentials together make the gravitational cavity for the atom (solid line), 
(right) We display the space-time quantum mechanical probability distribu- 
tion of the atomic wave packet in the Fermi accelerator represented as a 
quantum carpet. 

The study of the spatio-temporal dynamics of the atom in the gravita- 
tional cavity explains interesting dynamical features. In the long time dy- 
namics the atom undergoes self interference and exhibits revivals, and frac- 
tional revivals. Furthermore, as a function of space and time, the quantum 
interference manifests itself interestingly in the quantum carpets (GroBmann 
1997, Marzoli 1998), as shown in figure 4 (right). 

We show the quantum carpets for an atom in the Fermi accelerator, in 
figure 4. The dark gray regions express the larger probabilities, whereas, the 
in between light gray regions indicate smaller probabilities to find the atom in 
its space-time evolution. We recognize these regular structures not related to 
the classical space-time trajectories. Close to the surface of the atomic mirror, 
at z = 0, these structures appear as vertical canals sandwiched between two 
high probability dark gray regions. These canals become curved gradually 
away from the atomic mirror where the gravitational field is significant. 
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IV. 3. 2 Mode structure 

An atom observes almost an instantaneous impact with the atomic mirror, 
as it obeys two conditions: (i) The atom is initially placed away from the 
atomic mirror in the gravitational field, and (ii) the atomic mirror is made 
up of an exponentially decaying optical field with a very short decay length. 
Thus, we may take the gravitational cavity in a triangular well potential, 
made up of a linear gravitational potential and an infinite potential , where 
the atom undergoes a bounded motion 

We may express the corresponding effective Hamiltonian as 

H^^ + V(.), (26) 



where 



n^)-!""" '-n' (27) 



The solution to the stationary Schrodinger equation, Hipn = Enipn, hence 
provides 

ipn = J^ Ai{z - Zn) (28) 

as the eigen function (Wallis 1992). Here, A/" expresses the normalization 
constant. Furthermore, we take x = {2m'^g/h'^Y^^z as the dimensionless 
position variables, and Zn = {2m^g/h'^Y^^Zn as the nth zero of the Airy 
function. The index n therefore defines the nth mode of the cavity. Here 
the, nth energy eigen value is expressed as 

En = {-mh^g^/hn = mgzn. (29) 
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Figure 5: We express the Wigner distribution function for the first three 
eigen functions, defined in Eq. (28). 

In figure 5, we show the Wigner distribution, 
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W{z,p) = —- rn{z + y/2)M^~y/2)e'^ydy, (30) 

Inn j-oo 

for the first three eigen-functions, i.e. n = 1,2,3, in the scaled coordinates, 
z, and p. The distributions are symmetric around p = axes. Moreover, 
we note that the distributions extends in space along, 2;-axis as n increases. 
We can easily identify the non-positive regions of the Wigner distribution 
functions, as well. 

IV. 4 Optical traps 

A slight change of the shape of the atomic mirror from fiat to concave helps 
to make successful trap for atoms in a gravitational cavity (Wallis 1992, 
Ovchinnikov 1995). In addition, atomic confinement is possible by designing 
bi-dimensional light traps on a dielectric surface (Desbiolles 1996) or around 
an optical fiber (Fam Le Kien 2005). 

By appropriate choice of attractive and/or repulsive evanescent waves, we 
can successfully trap or guide atoms in a particular system. In the presence 
of a blue detuned optical field and another red detuned optical field with 
a smaller decay length, a net potential is formed on or around a dielectric 
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Figure 6: (left) Schematic diagram of a bi-dimensional atomic trap around 
an optical fiber, (right) Transverse plane profile of the total potential, Utot, 
produced by net optical potential and the van der Waals potential. (Fam Le 
Kien et al. 2004). 

surface. The bi-dimensional trap can be used to store atoms, as shown in 
figure 6. 

We may make an optical cylinder within a hollow optical fiber to trap 
and guide cold atoms (Renn 1995, Ito 1996). A laser light propagating in the 
glass makes an evanescent wave within the fiber. The optical field is detuned 
to the blue of atomic resonance, thus it exerts repulsive force on the atoms, 
leading to their trapping at the center of the fiber. The system also serves 
as a useful waveguide for the atoms. 

V Complex systems in atom optics 

In atom optics systems, discussed in section IV, the atomic dynamics takes 
place in two dimensions. However in the presence of an external driving force 
on the atoms an explicit time dependence appears in these systems. The 
situation may arise due to a phase modulation and/or amplitude modulation 
of the optical field. We may make the general Hamiltonian description of 
these driven systems as, 

H = Ho{x,p)-V{t)u{x + ^(t)). (31) 
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Here, Hq controls atomic dynamics in the absence of driving force and along 
X-axis. Furthermore, V{t) and ip{t) are periodic functions of time, and u 
mentions the functional dependence of the potential in space. 

Equation (31) reveals that the presence of explicit time dependence in 
these systems provides a three dimensional phase space. Hence, in the pres- 
ence of coupling, these systems fulfill the minimum criteria to expect chaos. 
Various situations have been investigated and important understandings have 
been made regarding the atomic evolution in such complex systems. Follow- 
ing we make a review of the atom optics systems studied in this regard. 

V.l Phase modulated standing wave field 

In 1992, Graham, Schlautmann and Zoller (GSZ) studied the dynamics of an 
atom in a monochromatic electromagnetic standing wave field made up of two 
identical and aligned counter propagating waves. As one of the running waves 
passes through an electro-optic modulator, a phase modulation is introduced 
in the field (Graham 1992). 

The effective Hamiltonian which controls the dynamics of the atom is, 

2 

H=^-Vocosikx + ipit)), (32) 

where, Vq expresses the constant amplitude and ip = xosm{Qt) defines the 
phase modulation of the field, with a frequency Q and amplitude xq. 

According to GSZ model, the atom in the phase modulated standing wave 
field experiences random kicks. Thus, a momentum is transferred from the 
modulated field to the atom along the direction of the field. 
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In classical domain the atom exhibits classical chaos as a function of the 
strength of the phase modulation and undergoes diffusive dynamics. How- 
ever, in the corresponding quantum dynamics the diffusion is sharply sup- 
pressed and the atom displays an exponentially localized distribution in the 
momentum space (Moore 1994, Bardroff 1995). 

V.2 Amplitude modulated standing wave field 

The atomic dynamics alters as the atom moves in an amplitude modulated 
standing wave field instead of a phase modulated standing wave field. The 
amplitude modulation may be introduced by providing an intensity modu- 
lation to the electromagnetic standing wave field through an acousto-optic 
modulator. 

The effective Hamiltonian which controls the atomic dynamics can now 
be expressed as 

2 

H = — -Vo cos(fit) cos(kx). (33) 

2m 

Thus, the system displays a double resonance structure as 

H = ^-Y^ [cos{kx + nt) + cos{kx - Qt)] . (34) 

The atom, hence, finds two primary resonances at +Q and —Q, where it 
rotates clockwise or counter clockwise with the field (Averbuckh 1995, Gorin 
1997, Monteoliva 1998). 
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V.3 Kicked rotor model 

The group of Mark Raizen at Austin, Texas presented the experimental re- 
ahzation of the Delta Kicked Rotor in atom optic (Robinson 1995, Raizen 
1999). This simple system is considered as a paradigm of chaos (Haake 
2001). In their experiment a cloud of ultra cold sodium atoms experiences 
a one dimensional standing wave field which is switched on instantaneously, 
and periodically after a certain period of time. 

The standing light field makes a periodic potential for the atoms. The 
field frequency is tuned away from the atomic transition frequency. There- 
fore, we may ignore change in the probability amplitude of any excited state 
as a function of space and time in adiabatic approximation. The general 
Hamiltonian which effectively governs the atomic dynamics in the ground 
state therefore becomes, 

H = Vo cos kx. (35) 

2m 

Here, the amplitude Vq is directly proportional to the intensity of the elec- 
tromagnetic field and inversely to the detuning. 

The simple one dimensional system may become non-integrable and dis- 
play chaos as the amplitude of the spatially periodic potential varies in time. 
The temporal variations are introduced as a train of pulses, each of a cer- 
tain finite width and appearing after a definite time interval, T. Thus, the 
complete Hamiltonian of the driven system appears as, 

2 +00 

H = ^—-Vocoskx J2 S{t-nT). (36) 



n=—oo 
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The atom, therefore, experiences a potential which displays spatial as well 
as temporal periodicity. 

It is interesting to note that between two consecutive pulses the atom 
undergoes free evolution, and at the onset of a pulse it gets a kick which 
randomly changes its momentum. The particular system, thus, provides an 
atom optics realization of Delta Kicked Rotor and enables us to study the 
theoretical predictions in a quantitative manner in laboratory experiments . 

V.4 Triangular billiard 

When a laser field makes a triangle in the gravitational field the atoms may 
find themselves in a triangular billiard. An interesting aspect of chaos enters 
depending on the angle between the two sides forming the billiard. Indeed 
the atomic dynamics is ergodic if the angle between the laser fields is an 
irrational multiple of vr (Artuso 1997a, Artuso 1997b) and may be pseudo 
integrable if the angles are rational (Richens 1981). 

V.5 An Atom optics Fermi accelerator 

The work horse of the Fermi accelerator is the gravitational cavity. In the 
atomic Fermi accelerator, an atom moves under the influence of gravitational 
field towards an atomic mirror made up of an evanescent wave field. The 
atomic mirror is provided a spatial modulation by means of an acousto-optic 
modulator which provides intensity modulation to the incident laser light 
field (Saif 1998). 

37 



Hence, the ultra cold two-level atom, after a normal incidence with the 
modulated atomic mirror, bounces off and travels in the gravitational field, 
as shown in figure 7. In order to avoid any atomic momentum along the 
plane of the mirror the laser light which undergoes total internal reflection, 
is reflected back. Therefore, we find a standing wave in the plane of the 
mirror which avoids any specular reflection (Wallis 1995). 

The periodic modulation in the intensity of the evanescent wave optical 
field may lead to the spatial modulation of the atomic mirror as 

I{z,t) = Ioexp[~2KZ + esm{nt)]. (37) 

Thus, the center-of-mass motion of the atom in z- direction follows effectively 
from the Hamiltonian 

H = ^-^ + mgz + hVleff exp[-2/tz + e sin(fit)], (38) 

where, f^e// denotes the effective Rabi frequency. Moreover, e and Vt express 
the amplitude and the frequency of the external modulation, respectively. 
In the absence of the modulation the effective Hamiltonian, given in equa- 
tion (38), reduces to equation (25). 

In order to simplify the calculations, we may make the variables dimen- 
sionless by introducing the scaling, z = {Q'^/g)z, p = {Q/mg)pz and t = Qt. 
Thus, the Hamiltonian becomes, 

2 

H = ^ + z + Voexp[-K{z- Xsint)] , (39) 

where, we take H = {fi'^/mg'^)H, Vq = TiVtef fVt^ I mg"^ , k = ng/Q'^, and 
A = eVtyg. 
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Figure 7: A cloud of atoms is trapped and cooled in a magneto-optic trap 
(MOT) to a few micro-Kelvin. The MOT is placed at a certain height above 
a dielectric slab. An evanescent wave created by the total internal reflec- 
tion of the incident laser beam from the surface of the dielectric serves as a 
mirror for the atoms. At the onset of the experiment the MOT is switched 
off and the atoms move under gravity towards the exponentially decaying 
evanescent wave field. Gravity and the evanescent wave field form a cavity 
for the atomic de Broglie waves. The atoms undergo bounded motion in this 
gravitational cavity. An acousto-optic modulator provides spatial modula- 
tion of the evanescent wave field. This setup serves as a realization of an 
atom optics Fermi accelerator. 

VI Classical chaos in Fermi accelerator 

The understanding of the classical dynamics of the Fermi accelerator is de- 
veloped together with the subject of classical chaos (Lieberman 1972, Licht- 
enberg 1980, 1983, 1992). The dynamics changes from integrable to chaotic 
and to accelerated regimes in the accelerator as a function of the strength of 
modulation. Thus, a particle in the Fermi accelerator exhibits a rich dynam- 
ical behavior. 

In the following, we present a study of the basic characteristics of the 
Fermi accelerator. 

VI. 1 Time evolution 

The classical dynamics of a single particle bouncing in the Fermi accelerator 
is governed by the Hamilton's equations of motion. The Hamiltonian of the 
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system, expressed in Eq. (39), leads to the equations of motion as 

.^-^p. (40) 

t> = -—- = -1 + KKoe-"^''^""*) . (41) 

oz 

In the absence of external modulation, the equations (40) and (41) are 
nonlinear and the motion is regular, with no chaos (Langhoff 1971, Gibbs 
1975, Desko 1983, Goodins 1991, Whineray 1992, Seifert 1994, Bordo 1997, 
Andrews 1998, Gea-Banacloche 1999). However, the presence of an exter- 
nal modulation introduces explicit time dependence, and makes the system 
suitable for the study of chaos. 

In order to investigate the classical evolution of a statistical ensemble in 
the system, we solve the Liouville equation (Gutzwiller 1992). The ensemble 
comprises a set of particles, each defined by an initial condition in phase 
space. Interestingly, each initial condition describes a possible state of the 
system. 

We may represent the ensemble by a distribution function, Pq{zo,pq) at 
time t = in position and momentum space. The distribution varies with 
the change in time, t. Hence, at any later time t the distribution function 
becomes, P{z,p,t). 

In a conservative system the classical dynamics obeys the condition of 
incompressibility of the flow (Lichtenberg 1983), and leads to the Liouville 
equation, expressed as, 

40 



Here, p is the force experienced by every particle of the ensemble in the Fermi 
accelerator, defined in equation (41). Hence, the classical Liouville equation 
for an ensemble of particles becomes, 

— +p— - [1 - KVoe^Y>{-<z - Asint)}]— \ P{z,p,t) = . (43) 

The general solution of equation (43) satisfies the initial condition, P{z,p, t 
0) = Pq{zq,po). According to the method of characteristics (Kamke 1979), it 
is expressed as 

/OO /'OO 

dzo / dpo6{z - z{zQ,po,t)}6{p - p{zo,po,t)}Po{zo,Po) ■ 
-OO J —OO 

(44) 
Here, the classical trajectories z = z{zQ,pQ,t) and p = p{zQ,pQ,t) are the 
solutions of the Hamilton equations of motion, given in equations (40) and 
(41). This amounts to say that each particle from the initial ensemble follows 
the classical trajectory {z,p). As the system is nonintegrable in the presence 
of external modulation, we solve equation (44) numerically. 

VI. 2 Standard mapping 

We may express the classical dynamics of a particle in the Fermi accelera- 
tor by means of a mapping. The mapping connects the momentum of the 
bouncing particle and its phase just before a bounce to the momentum and 
phase just before the previous bounce. This way the continuous dynamics of 
a particle in the Fermi accelerator is expressed as discrete time dynamics. 
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In order to write the mapping, we consider that the modulating surface 
undergoes periodic oscillations following sinusoidal law. Hence, the position 
of the surface at any time is z = Asint, where A defines the modulation 
amplitude. In the scaled units the time of impact, t, is equivalent to the 
phase ip. 

Furthermore, we consider that the energy and the momentum remain 
conserved before and after a bounce and the impact is elastic. As a result, 
the bouncing particle gains twice the momentum of the modulated surface, 
that is, 2Acos(/9, at the impact. Here, we consider the momentum of the 
bouncing particle much smaller than that of the oscillating surface. Moreover, 
it undergoes an instantaneous bounce. 

Keeping these considerations in view, we express the momentum, Pi+i, 
and the phase, v^i+i, just before the {i + l)th bounce in terms of momentum. 
Pi, and the phase, ipi, just before the ith bounce, as 

Pi+i = -pi - Ati + 2A cos (fi , 

^i+i = ipi + Ati • (45) 

Here, the time interval, Ati = ^i+i ~ ^i, defines the time of flight between 
two consecutive bounces. Hence, the knowledge of the momentum and the 
phase at the ith bounce leads us to the time interval Ati as the roots of the 
equation, 

PiAti - -Atf = A(sin(v9i + Ati) - siny^i) . (46) 
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We consider that the amphtude of the bouncing particle is large enough 
compared to the amplitude of the external modulation, therefore, we find no 
kick-to-kick correlation. Thus, we may assume that the momentum of the 
particle just before a bounce is equal to its momentum just after the previous 
bounce. The assumptions permit us to take, Ati ~ — 2pi+i. As a result the 
mapping reads, 

Pi+i =Pi-2Acos(/9i , 

V?i+i = v^i - 2pi+i . (47) 

We redefine the momentum as pi = —2p[, and K = 4A. The substitutions 
translate the mapping to the standard Chirikov- Taylor mapping (Chirikov 
1979), that is, 

pi+i = pi + K cos ipi , 

V^i+i = Vi + Pi+i ■ (48) 

Hence, we can consider the Fermi accelerator as a discrete dynamical system 
(Kapitaniak 2000). 

The advantage of the mapping is that it depends only on the kick strength 
or chaos parameter, K = 4A. Therefore, simply by changing the value of the 
parameter, K, the dynamical system changes from, stable with bounded 
motion to chaotic with unbounded or diffusive dynamics. There occurs a 
critical value of the chaos parameter, at which the change in the dynamical 
characteristics takes place, is K = Kcr = 0.9716.... (Chirikov 1979, Greene 
1979). 
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VI. 3 Resonance overlap 

Periodically driven systems, expressed by the general Hamiltonian given in 
equation (31), exhibit resonances. These resonances appear whenever the fre- 
quency of the external modulation, Q, matches with the natural frequency 
of unmodulated system, uj (Lichtenberg 1992, Reichl 1992). Thus, the reso- 
nance condition becomes, 

Niu- Mn = 0, (49) 

where, A^ and M are relative prime numbers. These resonances are spread 
over the phase space of the dynamical systems. 

Chirikov proved numerically that in a discrete dynamical system ex- 
pressed by Chirikov mapping, the resonances remain isolated so far as the 
chaos parameter, K, is less than a critical value, Kcr = 1 (Chirikov 1979). 
Later, by his numerical analysis, Greene established a more accurate measure 
for the critical chaos parameter as K^r — 0.9716... (Greene 1979). Hence, the 
dynamics of a particle in the system remains bounded in the phase space by 
Kolmogorov-Arnold-Moser (KAM) surfaces (Arnol'd 1988, Arnol'd 1968). As 
a result, only local diffusion takes place. 

Following our discussion presented in the section VI. 2, we note that in the 
Fremi accelerator the critical chaos parameter. A/, is defined as A^ = -ft'cr/4 — 
0.24. Hence, for a modulated amplitude much smaller than A; the phase space 
is dominated by the invariant tori, defining KAM surfaces. These surfaces 
separate resonances. However, as the strength of modulation increases, area 
of the resonances grow thereby more and more KAM surfaces break. 
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At the critical modulation strength, A;, last KAM surface corresponding 
to golden mean is broken (Lichtenberg 1983). Hence, above the critical mod- 
ulation strength the driven system has no invariant tori and the dynamics of 
the bouncing particle is no more restricted, which leads to a global diffusion. 

The critical modulation strength A/, therefore, defines an approximate 
boundary for the onset of the classical chaos in the Fermi accelerator. The dy- 
namical system exhibits bounded motion for a modulation strength A smaller 
than the critical value A;, and a global diffusion beyond the critical value. 

VI. 4 Brownian motion 

As discussed in section VI.2, classical dynamics of a particle in the Fermi 
accelerator is expressed by the standard mapping. This allows to write the 
momentum as 

3 

Pj = Po + 4A ^ cos ipn, (50) 

ra=l 

at the jth bounce. Here, po is the initial scaled momentum. In order to cal- 
culate the dispersion in the momentum space with the change of modulation 
strength. A, we average over the phase, ip. This yields, 

j j 
Ap/ = (p/) - {Pjf = (4A)^ £ £ (cos(v9„) cos(^„0) • (51) 

n=l n'=l 

As the amplitude of the bouncing particle becomes very large as com- 
pared with that of the oscillating surface, we consider that no kick-to-kick 
correlation takes place. This consideration leads to a random phase for the 
bouncing particle in the interval [0, 2tt] at each bounce above the critical 
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modulation strength, A/. Thus a uniform distribution of the phase appears 
over many bounces. Therefore, we get, 

(C0S(V9„) C0S(V9„/)) = {COS^ If n)5n,n' = -^K,n> ■ (52) 

Hence, equation (51) provides 

Ap^ = 2jX^=jD, (53) 

where D = 2A^ and expresses the diffusion constant. Moreover, j describes 
the number of bounces. 

Since we can find an average period over j bounces, the number of bounces 
grows hnearly with time. This leads to a linear growth of the square of width 
in momentum space as a function of the evolution time, as we find in the 
Brownian motion. 

It is conjectured that in long time limit the diffusive dynamics of the Fermi 
accelerator attains Boltzmann distribution. We may express the distribution 
as, 

PM,P) = (27r)-i/V'/'exp[-(pV2 + ^)/r7] , (54) 

where, the quantity rj represents effective temperature (Saif 1998). 

Hence, the momentum distribution follows Maxwellian distribution in the 
classical phase space, and is Gaussian, that is, 

P4p) = -l=exp[-py{2r])]. (55) 
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Moreover, the classical position distribution follows exponential barometric 
formula 

Pd{z) = -exp{-z/r]) , (56) 

predicted from Eq. (54). The conjecture is confirmed from the numerical 
study, as we shall see in section VIII. 1. The equations (55) and (56) lead to 
the conclusion that Ap^ = Az = jD. 

VI. 5 Area preservation 

Preservation of area in phase space is an important property of the Hamilto- 
nian systems (Pustyl'nikov 1977, Lichtenberg 1983) and corresponds to the 
conservation of energy. In a system the determinant of Jacobian, defined as, 

in, i^] ^ (57) 

leads to the verification of area preservation in a classical system. 

In the case of a non-dissipative, conservative system the determinant of 
the Jacobian is equal to one which proves the area-preserving nature of the 
system in phase space. However, for a dissipative system it is less than one. 

For the Fermi accelerator the determinant of the Jacobian, obtained with 
the help of equation (48), is 

J = det(^-^''^^\^]=l. (58) 

\ —K sm (/9i 1 / ^ ' 

Hence, the dynamics of a particle in the Fermi accelerator exhibits the prop- 
erty of area preservation in phase space. 
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VI. 6 Lyapunov exponent 

The exponential sensitivity of a dynamical system on the initial conditions 
is successfully determined by means of Lyapunov exponent. The measure of 
Lyapunov exponent is considered to determine the extent of classical chaos 
in a chaotic system (Schuster 1989, Gutzwiller 1992, Ott 1993). 
We may define the exponent as. 

Here, d{0) represents distance between two phase points at time, t = 0, and 
d{t) describes their distance after an evolution time, t. If the dynamics is 
regular and non diffusive, the exponent will be zero. However if it is diffusive 
we find nonzero positive Lyapunov exponent. A system displaying dissipative 
dynamics has a nonzero negative exponent. 

For weak modulation, that is for A < A^, we find the exponent zero in gen- 
eral. This indicates a dominantly regular dynamics in the Fermi accelerator 
for weak modulation strength. For modulation strengths A > A;, the dy- 
namics is dominantly diffusive in the system established by nonzero positive 
exponent (Saif 1998). 

VI. 7 Accelerating modes 

The classical work of Pustyl'nikov on the Fermi accelerator model (Pustyl'nikov 
1977) guarantees the existence of a set of initial data, such that, trajectories 
which originate from the set always speed-up to infinity. The initial data 
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appears always in circles of radius p > and has positive Lebesgue measure. 
Furthermore, the existence of these accelerating modes requires the presence 
of certain windows of modulation strength, A. 

The windows of modulation strength which supports accelerated trajec- 
tories, read as. 



STT < A < 1/1 + (s7r)2 . (60) 

Here, s takes integer and half integer values for the sinusoidal modulation of 
the reflecting surface. 

In order to probe the windows of modulation strength which support 
the unbounded acceleration, we may calculate the width of momentum dis- 



tribution, Ap = \l {p^) — {p) , numerically as a function of the modulation 
strength. Here, {p) and {p^) are the first and the second moment of momen- 
tum, respectively. 

We consider an ensemble of particles which is initially distributed fol- 
lowing Gaussian distribution. In order to study the dynamical behavior of 
the ensemble we record its width in the momentum space after a propaga- 
tion time for different modulation strengths. We note that for very small 
modulation strengths, the widths remain small and almost constant, which 
indicates no diffusive dynamics, as shown in figure 8. For larger values of A 
the width Ap increases linearly, which follows from the equation (53). 

We find that at the modulation strengths which correspond to the win- 
dows on modulation strengths, given in equation (60), the diffusion of the 
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Figure 8: The width of the momentum distribution Ap is plotted as a func- 
tion of the modulation strength, A. An ensemble of particles, initially in 
a Gaussian distribution, is propagated for a time, t = 500. The numeri- 
cal results depict the presence of accelerated dynamics for the modulation 
strengths, expressed in equation (60). 

ensemble is at its maximum. The behavior occurs as the trajectories which 
correspond to the areas of phase space supporting accelerated dynamics un- 
dergo coherent acceleration, whereas the rest of the trajectories of the initial 
distribution display maximum diffusion. For the reason, as the modulation 
strengths increase beyond these values, the dispersion reduces as we find in 
Fig. 8. 

Equation (60) leads to the modulation strength. Am, for which maximum 
accelerated trajectories occur in the system. We find these values of A^ as. 



\m = ^, . (61) 



The value is confirmed by the numerical results (Saif 1999), as expressed in 
figure 8. 

For a modulation strength given in equation (60) and an initial ensemble 
originating from the area of phase space which supports accelerated dynam- 
ics, we find sharply suppressed value of the dispersion. It is a consequence of 
a coherent acceleration of the entire distribution (Yannocopoulos 1993, Saif 
1999). 
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VII Quantum dynamics of the Fermi acceler- 
ator 

The discrete time dynamics of a classical particle in the Fermi accelerator 
is successfully described by the standard mapping. For the reason, we find 
onset of chaos in the system as the chaos parameter exceeds its critical value. 
The presence of classical chaos in the Fermi accelerator is established by non- 
zero positive Lyapunov exponent (Saif 1998). This makes the dynamics of 
a quantum particle in the Fermi accelerator quite fascinating in the frame 
work of quantum chaos. 

Similar to classical dynamics, we find various dynamical regimes in the 
quantum evolution of the system as a function of modulation strength. The 
quantum system however has another controlling parameter, that is, the 
scaled Planck's constant. The commutation relation between the dimension- 
less coordinates, z and p, defined for equation (39), appear as 

[z,p]=i^. (62) 

This leads to the definition of the scaled Planck's constant as fe = (r2/r2o)^. 
The frequency f^o takes the value as f^o = {mg^ ITif-l"^ . Thus, it is easily 
possible to move from the semi-classical to pure quantum mechanical regime 
simply by changing the frequency of the external modulation alone. 

The quantum dynamics of the center-of-mass motion of an atom in the 
Fermi accelerator, follows from the time dependent Schrodinger equation. 



'^-Qt 



p2 

h -2 + Vo exp [— /t(-2 — A sint)] 



^ . (63) 
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The equation (63) describes the dynamics of a cold atom which moves to- 
wards the modulated atomic mirror under the influence of gravitational field 
and bounces off. In its long time evolution the atom experiences large number 
of bounces. 

A quantum particle, due to its nonlocality, always experiences the mod- 
ulation of the atomic mirror during its evolution in the Fermi accelerator. 
However, the classical counterpart "feels" the modulation of the mirror only 
when it bounces off. This contrast in quantum and classical evolution leads 
to profound variations in the dynamical properties in the two cases. In the 
next sections we study the quantum characteristics of the Fermi accelerator. 

VII. 1 Near integrable dynamics 

As we discussed in section VI. 3, for a modulation strength A smaller than A^, 
classical resonances remain isolated. The classical dynamics of a particle in 
the Fermi accelerator is, therefore, restricted by KAM surfaces. As a result, 
after its initial spread over the area of resonances, a classical ensemble stops 
diffusing and the classical dynamics stays bounded. 

Similar to the isolated classical resonances, isolated quantum resonances 
prone to exist for smaller modulations (Berman and Kolovsky 1983, Reichel 
and Lin 1986). For the reason an evolution similar to the classical evolution is 
found even in the quantum domain. We find this behavior as a wave packet in 
its quantum dynamics mimics classical bounded motion, both in the position 
and momentum space in the Fermi accelerator, for smaller modulations. 
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Manifestation of the behavior comes from the saturation of the width, in 
the position space I\z and momentum space Ap, as a function of time, both 
in the classical and quantum domain. The saturation values for the widths 
in the position space and the momentum space may differ due to the size of 
underlying resonances (Saif et al. 1998, Saif 1999). 

VII. 2 Localization window 

In the Fermi accelerator the classical system undergoes a global diffusion 
as the resonance overlap takes place above the critical modulation strength, 
A; = 0.24. In contrast, in the corresponding quantum mechanical domain 
there occurs another critical modulation strength, A„, which depends purely 
on the quantum laws. 

At the critical modulation strength, A„, a phase transition occurs and the 
quasi-energy spectrum of the Floquet operator changes from a point spectrum 
to a continuum spectrum (Benvenuto 1991, Oliveira 1994, Brenner 1996). 
We may define the critical modulation strength for the Fermi accelerator as, 
A„ = vfc/2, when the exponential potential of the atomic mirror is considered 
as an infinite potential (Saif 1998). Beyond this critical modulation strength 
quantum diffusion takes place. 

Hence, the two conditions together make a window on the modulation 
strength. We may find a drastic difference between the classical dynamics 
and the corresponding quantum dynamics, for a modulation strength which 
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fulfills the condition, 

A^ < A < A„ . (64) 

Within the window classical diffusion sets in whereas the corresponding 
quantum dynamics displays localization. For the reason, we name it as lo- 
calization window (Chen f997, Saif f998, Saif f999). 

VII. 3 Beyond localization window 

Above the upper boundary of the localization window, A„, even the quan- 
tum distributions display diffusion. For the reason, a quantum wave packet 
maintains its widths, in the momentum and the position space, only up to 
A ~ A„. Beyond the value it starts spreading similar to the classical case. 

For the reason, above the critical modulation strength, X^, width in the 
momentum space Ap and in the position space Az, display a growing be- 
havior. However, in addition to an overall growing behavior, there exist 
maxima and minima in a regular fashion similar to the classical case, as the 
modulation strength A increases. We show the classical behavior in figure 8. 

The presence of maxima corresponds to maximum dispersion. Interest- 
ingly, the maxima occur at A = Am as expressed in equation (61), and the 
size of the peaks is determined from equation (60). Hence, we infer that the 
behavior is a consequence of the accelerated trajectories, and in accordance 
with the Pustyl'nikov's work, as discussed in section VI. 7. 

For a modulation strength given in equation (60), the portions of a prob- 
ability distribution which originate from the area of phase space supporting 
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Figure 9: The classical and the quantum dynamics are compared far above 
the localization window by calculating the momentum distributions, P{p), 
after the evolution time, t = 500. (left) For A = 1.7 and (right) for A = 2.4, 
marked in Fig. 8, we plot the momentum distributions in both the classical 
and the quantum space, together, as mirror images. We find spikes appearing 
in the momentum distributions for A = 1.7, which is due to the presence of 
coherent accelerated dynamics. 

accelerated dynamics, always get coherent acceleration. However, the rest 
of the distribution displays maximum diffusion and thus attains maximum 
widths. 

The coherent acceleration manifests itself in the regular spikes in the mo- 
mentum and the position distributions, as shown in figure 9. It is important 
to note that the spiky behavior takes place for the modulation strengths 
which satisfy the condition given in equation (60). 

These spikes gradually disappear as we choose the modulation strength A 
away from these windows. The behavior is a beautiful manifestation of the 
quantum non- dispersive accelerated dynamics in the Fermi accelerator. 

VIII Dynamical localization of atoms 

The beautiful work of Casati and co-workers on the delta kicked rotor, now a 
paradigm of classical and quantum chaos, led to the discovery of dynamical 
localization (Casati 1979). They predicted that a quantum particle follows 
diffusive dynamics of a classical chaotic system only up to a certain time, 
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named as quantum break time. Beyond the time the classical diffusion is 
ceased due to the quantum interference and the initial quantum distribution 
settles into an exponential localization. 

Similar behavior is predicted to occur in the quantum dynamics of an 
atom in the Fermi accelerator (Chen 1997, Saif 1998). However, the exis- 
tence of the dynamical localization is restricted to the localization window 
(Benvenuto 1991), defined in the equation (64). Furthermore, the effective 
Planck's constant imposes another condition as it is to be larger than a min- 
imum value, 4Af. 

VIII. 1 Probability distributions: An analysis 

An atomic wave packet, expressed initially as the ground state of harmonic 
oscillator and satisfying minimum uncertainty criteria, displays an exponen- 
tial localization behavior when propagated in the Fermi accelerator. The 
localization takes place both in the momentum space and position space. 
Furthermore, it occurs for a time beyond quantum break time, and for a 
modulation within the localization window. The final quantum distributions 
are a manifestation of interference phenomena, thereby independent of the 
choice of initial distributions of the wave packet. 

In the momentum space the atomic wave packet redistributes itself around 
the initial mean momentum after an evolution time of many bounces (Lin 
1988). In the long time limit, the probability distribution decays exponen- 
tially in the momentum space, as shown in figure 10. It is estimated that the 
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probability distribution in momentum space follows, exp (— |p|/£), behavior. 
Here, d, describes the localization length. 

In the corresponding position space the wave packet reshapes itself in 
a different manner. The mean position of the probability distribution is 
estimated to be still at the initial mean position of the wave packet. However, 
it decays exponentially in the gravitational field when the exponent follows 
the square root law (Saif 2000a), as shown in the figure 10. 

Thence, the quantum distributions in the position space and in the mo- 
mentum space are completely different from their classical counterparts, dis- 
cussed in section 6.4. In contrast to the quantum mechanical exponential 
localization in momentum space, classical ensemble undergoes diffusion and 
attains Maxwell distribution, defined in equation (55). 

Moreover in the position space the quantum distribution is exponential 
with the exponent following square root behavior, whereas, it follows expo- 
nential barometric equation in the classical domain and the exponent follows 
linear behavior, as given in equation (56). Thus, the quantum distribution 
displays much rapid decay in the position space. 

VIII. 2 Dispersion: Classical and quantum 

In the quantum dynamics of the delta kicked rotor model, energy grows 
following the classical diffusive growth only up to quantum break time, es- 
timated as t* ^ A^/fc . Later, it saturates and deviates from the unlimited 
growth in energy which takes place in the classical case. 
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Figure 10: Comparison between the classical (thick lines) and the quantum 
mechanical (thin lines) distributions in momentum space (left) and in posi- 
tion space (right) at a fixed evolution time: We have k = 4 for which the 
localization window reads 0.24 < A < 1. We take the modulation strength, 
A = 0.8, which is within the window. The initial minimum uncertainty wave 
packet attains exponential distribution after a scaled evolution time t = 3200 
in the Fermi accelerator, with the exponents following linear and square root 
behavior in momentum space and in position space, respectively. The cor- 
responding classical distributions follow Maxwell distribution and baromet- 
ric equation in momentum space and in position space, respectively. Here, 
Az = 2 and the number of particles in the classical simulation is 5000. The 
time of evolution, t = 3200, is far above the quantum break time t* = 250. 
The estimated behavior is represented by thin dashed lines. 

In a dynamical system the suppression of classical diffusion in quantum 
domain is regarded as an evidence of dynamical localization (Reichl 1986). 
It is conjectured as independent of initial distribution (Izrailev 1990). 

G. J. Milburn discussed that a wave packet displays saturation in its en- 
ergy after quantum break time in the Fermi accelerator, provided the mod- 
ulation strength satisfies the condition given in equation (64) (Chen and 
Milburn 1997). This saturation in energy, when the corresponding classical 
system supports unlimited gain, establishes the quantum suppression of the 
classical chaos in the Fermi accelerator, as shown in figure 11. 

In the Fermi accelerator the phenomenon of dynamical localization exists 
in the position and momentum space, as well. In the classical dynamics of an 
ensemble in the Fermi accelerator, the dispersion in the position space and the 
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Figure 11: Comparison between classical (dashed line) and quantum (solid 
line) dynamics within the localization window: The classical and quantum 
evolution go together till the quantum break time, t*. Above the quantum 
break time, while the classical ensemble diffuses as a function of time, the 
quantum system shows a suppression of the diffusion. (Chen 1997). 

square of the dispersion in the momentum space are calculated numerically. 
It is observed that they exhibit a linear growth as a function of time. 

It is straightforward to establish the behavior analytically as Ap^ ~ Az ~ 
Dt, with the help of equations (55) and (56). The corresponding quantum 
mechanical quantities, however, saturate or display significantly slow increase 
(Saif 1998, Saif 2000b). 

VIII. 3 Effect of classical phase space on dynamical lo- 
calization 

As long as isolated resonances exist in the classical phase space, the localiza- 
tion phenomenon remains absent in the corresponding quantum dynamics. 
The dynamical localization however appears as the overlap of the resonance 
takes place in the classical dynamics. As a consequence in the quantum do- 
main, the probability distributions decay exponentially in the stochastic sea 
(Chirikov and Shepelyansky 1986). 

In the Fermi accelerator the quantum probability distributions display 
exponential localization. However, the distribution in position space exhibits 
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a different decay, as the exponent follows square root behavior in contrast to 
a linear behavior as in the momentum space (Saif et al. 1998, Saif 1999). 

We note that in addition to an overall decay there occur plateau struc- 
tures in the two distributions both in the classical and in the quantum dy- 
namics. The numerical results provide an understanding to the behavior as 
the confinement of the initial probability distributions. 

A portion of the initial probability distribution which originates initially 
from a resonance finds itself trapped there. For the reason over the area of a 
resonance the approximate distribution of the probability is uniform which 
leads to a plateau (Chirikov 1986, Saif 2000b). 

In the presence of the numerical study we conjecture that in the quantum 
domain corresponding to a classical resonance there occur a discrete quasi- 
energy spectrum, whereas in the stochastic region we find a quasi continuum 
of states (Saif 2000c). For the reason, the quantum probability distribution 
which occupies the quasi-continuum spectrum undergoes destructive inter- 
ference and manifests the exponential decay of the probability distributions. 

The height of plateaus appearing in the probability distributions rises 
for the higher values of the Planck's constant. However, their size and the 
location remain the same (Lin 1988). Interestingly, for the larger values of 
the effective Planck's constant, k, the number of plateaus increases as well. 

These characteristics are comprehensible as we note that, for a larger val- 
ues of k, the size of thee initial minimum uncertainty wave packet increases. 
For the reason, the portion of the initial probability distribution which enters 
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the stable islands increases, and therefore, displays an increase in the height 
of the plateaus. Moreover, for a larger values of k, the exponential tails of 
the quantum distributions cover more resonances leading to more plateau 
structures in the localization arms (Saif 2000b). 

VIII. 4 Quantum delocalization 

Above the upper boundary of the localization window, the interference is 
overall destructive and leads to the quantum mechanical delocalization in 
the Fermi accelerator. Thus, we find exponential localization within the 
localization window only, as shown in figure 12. We may predict the presence 
of quantum delocalization in a time dependent system by observing the time 
independent spectrum of the potential. 

In nature there exist two extreme cases of potentials: Tightly binding po- 
tentials, for which the level spacing always increases between adjacent energy 
levels as we go up in the energy, for example, in one dimensional box; Weakly 
binding potentials, for which the level spacing always reduces between the 
adjacent levels, and we find energy levels always at a distance smaller than 
before, as we go up in the energy, for example, in the gravitational potential. 
Quadratic potential offers a special case, as the level spacing remains equal 
for all values of energy. Therefore, it lies between the two kinds of potentials. 

The fundamental difference of the time independent potentials has a dras- 
tic influence on the properties of the system in the presence of external pe- 
riodic driving force (Saif 1999). In the presence of an external modulation, 
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Figure 12: A transition from a localized to a delocalized quantum mechanical 
momentum distribution occurs at A^ = 1, for k = 4. (a) For the modulation 
strength A = 0.8 < A„, we find exponential localization, (b) In contrast for 
A = 1.2 > A„ we find a broad Gaussian distribution indicating delocalization, 
similar to Maxwell distribution in classical diffusion. The initial width of the 
atomic wave packet in the Fermi accelerator is, Az = 2 in position space, 
and Ap = k /2Az = 1 in momentum space. We take V = 4 and k, = 0.5. 
The quantum distributions are noted after a propagation time t = 3200. 

for example, the spectrum of a tightly binding potential is a point spectrum 
(Hawland 1979, 1987, 1989a, 1989b, Joye 1994). 

In the presence of an external modulation, we observe a transition from 
point spectrum to a continuum spectrum in the weakly binding potentials, 
above a certain critical modulation strength (Delyon 1985, Oliveira 1994, 
Brenner 1996). Therefore, in the weakly binding potentials we find local- 
ization below a critical modulation strength, where the system finds a point 
spectrum. 

The quantum delocalization occurs in the Fermi accelerator as a material 
wave packet undergoes quantum diffusion for a modulation strength above 
Xu- In this case the momentum distribution becomes Maxwellian as we find 
it in the corresponding classical case. The presence of the diffusive dynamics 
above the upper boundary, shown in figure 12(b), is an illustration of the 
quasi- continuum spectrum in the Fermi accelerator above A„. 
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IX Dynamical recurrences 

The quantum recurrence phenomena are a beautiful combination of classical 
mechanics, wave mechanics, and quantum laws. A wave packet evolves over 
a short period of time in a potential, initially following classical mechanics. 
It spreads while moving along its classical trajectory, however rebuilds itself 
after a time called as classical period. In its long time evolution it follows 
wave mechanics and gradually observes a collapse. However, the discreteness 
of quantum mechanics leads to the restoration and restructuring of the wave 
packet at the time named as quantum revival time. 

In one degree of freedom systems the phenomena of quantum revivals 
are well studied both theoretically and experimentally. The quantum re- 
vivals were first investigated in cavity quantum electrodynamics (Eberly 
1980, Narozhny 1981, Yurke 1986). Recently, the existence of revivals has 
been observed in atomic (Parker 1986, Alber 1986, Alber 1988, Averbukh 
1989, 1989a, Braun 1996, Leichtle 1996, 1996a, Aronstein 2000) and molec- 
ular (Fischer 1995, Vrakking 1996, Grebenshchikov 1997, Doncheski 2003) 
wave packet evolution. 

The periodically driven quantum systems (Hogg 1982, Haake 1992, Bres- 
lin 1997), and two-degree-of- freedom systems such as stadium billiard (Tomso- 
vic and Lefebvre 1997) presented the first indication of the presence of quan- 
tum revivals in higher dimensional systems. Later, it is proved that the 
quantum revivals are generic to the periodically driven one degree of freedom 
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quantum systems (Saif 2002, Saif 2005) and to the coupled two degrees-of- 
freedom quantum systems (Saif 2005a). 

IX. 1 Dynamical recurrences in a periodically driven 
system 

The time evolution of a material wave packet in a one dimensional system 
driven by an external periodic field is governed by the general Hamiltonian, 
H, expressed in the dimensionless form as 

H = Ho + \V{z)sint . (65) 

Here, Hq, is the Hamiltonian of the system in the absence of external driving 
field. Moreover, A expresses the dimensionless modulation strength. We 
may consider |n), and. En as the eigen states and the eigen values of Hq, 
respectively. 

In order to make the measurement of the times of recurrence in the clas- 
sical and the quantum domain, we solve the periodically driven system. In 
the explicitly time dependent periodically driven system the energy is no 
more a constant of motion. Therefore, the corresponding time dependent 
Schrodinger equation is solved by using the secular perturbation theory in 
the region of resonances (Born 1960). 

In this approach, faster frequencies are averaged out and the dynamical 
system is effectively reduced to one degree of freedom. The reduced Hamil- 
tonian is integrable and its eigen energies and eigen states are considered as 
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the quasi-energy and the quasi-energy eigen states of the periodically driven 

system. 

IX. 1.1 Quasi-energy and quasi-energy eigen states 

In order to solve the time dependent Schrodinger equation which controls the 
evolution in a periodically driven system, Berman and Zaslavsky (Berman 
1977), and later Flatte' and Holthaus (Flatte' 1996), suggested its solution 
in an ansatz, as 



\m) = J2Cnit)\n)exp 



Er + (n- r] 



h 

N 



h 



(66) 



Here, h is the scaled Planck's constant for the driven system, Er is the energy 
corresponding to the mean quantum number r, and A^ describes the number 
of resonance. 

By following the method of secular perturbation the time dependent 
Schrodinger equation reduces to 



ih 



dg_ 
dt 



h'NX d' 



- iNh u 



d 



+ Ho + \V sm{9) 



g{9,t). (67) 



2 (9^2 --y- ^JQQ 

The parameters, C, = E"/h , and, cu = E'JU^ express the nonlinearity in the 
system and the classical frequency, respectively. Here, E'^ defines the second 
derivative of energy with respect to quantum number n calculated at n = r, 
and E'^ defines the first derivative of the unperturbed energy, calculated at 
n = r. 
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The function g{6) is related to C„(t) as, 



ZTT J 

^(^)e-i(n-'-)e/Af ^0 _ (68) 
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We may define g{6,t) as 

g[e^ t) = ^(^)e-*^*/^e-^(--i/^)^/^^" (69) 

and take the angle variable ^, as ^ = 2z + 7i/2. These substitutions reduce the 
equation (67) to the standard Mathieu equation (Abramowitz 1992), which 
is, 



.2 



& 

7— r + a- 2qcos(2z) 



~g{z) = . (70) 



Here, 

and q = AW / N"^ C^Ti^ . Moreover, gj^{z) is a 7r-periodic Floquet function, and 
Ey defines the quasi-energy of the system for the index z/, given as 

_ 2(n-r) 2(a; - 1/iV) 

Hence, the quasi energy eigen values £y, are defined as 

^^ = — ^a^(g) ^T + ^0 • (73) 

In this way, we obtain an approximate solution for a nonlinear resonance of 
the explicitly time-dependent system. 
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IX. 1.2 The dynamical recurrence times 

An initial excitation produced, at n = r, observes various time scales at which 
it recurs completely or partially during its quantum evolution. In order to 
find these time scales at which the recurrence occur in a quantum mechanical 
driven system, we employ the quasi energy, S^, of the system (Saif 2002, Saif 
2005). 

These time scales T^ , are inversely proportional to the frequencies, uj^^\ 
such that T^^ = 2tx /uj^'-'\ where j is an integer. We may define the frequency, 
uj'^^\ of the reappearance of an initial excitation in the dynamical system as 

calculated at n = r. The index j describes the order of differentiation of the 
quasi energy, £k with respect to the quantum number, n. The equation (74) 
indicates that as value of j increases, we have smaller frequencies which lead 
to longer recurrence times. 

We substitute the quasi energy, E^ as given in equation (73), in the ex- 
pression for the frequencies defined in the equation (74). This leads to the 
classical period, T{ = T^ for the driven system as, 

T^^') = [1 - M^^')]To^''^A, (75) 

and, the quantum revival time T^ = T^ , as, 

Tf ) = [1 - M^'^^]tS^\ (76) 
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Here, the time scales Tq and Tq express the classical period and the 
quantum revival time in the absence of external modulation, respectively. 
They are defined as Tq\= 2-k/uj), and Tq \= 2-k{^Q^^). Furthermore, 

A = (1 — cJAf/cj)^-*-, where Un = l/N. 

The time modification factor Mjf''^ and MJ^^^ are given as, 

M(-0 = -lf^^^V ' , (77) 



and 



where. 



M(«) = if^^^V^±i^, (78) 



Equations (75) and (76) express the classical period and the quantum 
revival time in a periodically system. These time scales are function of the 
strength of the periodic modulation A and the matrix element, V. Moreover, 
they also depend on the frequency, u, and the nonlinearity, (, associated 
with the nndriven or unmodulated system. 

As the modulation term vanishes, that is A = 0, the modification terms 
Mjf^^ and MJ^^^ disappear. Since there are no resonances for A = 0, we find 
A = 1. Thus, from equations (75) and (76), we deduce that in this case 
the classical period and the quantum revival time in the presence and in the 
absence of external modulation are equal. 
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IX. 2 Classical period and quantum revival time: In- 
terdependence 

The nonlinearity ( present in the energy spectrum of the nndriven system, 
contributes to the classical period and the quantum revival time, in the 
presence and in the absence of an external modulation. Following we analyze 
different situations as the nonlinearity varies. 

IX. 2.1 Vanishing nonlinearity 

In the presence of a vanishingly small nonlinearity in the energy spectrum, 
i.e for C ~ 0, the time modification factors for the classical period Mjf'-^ and 
the quantum revival time M^*^^ vanish, as we find in equations (77) and (78). 
Hence, the modulated linear system displays recurrence only after a classical 
period, which is TJ'^ = Tg A = 27rA/c(j. The quantum revival takes place 
after an infinite time, that is, T^ = Tq = oo. 

IX. 2. 2 Weak nonlinearity 

For weakly nonlinear energy spectrum, the classical period, T^^ , and the 
quantum revival time, TJ , in the modulated system are related with, Tg , 
and, Tg , of the linmodulated system as. 

As discussed earlier, the quantum revival time, Tg , depends inversely 
on the nonlinearity, (, in the unmodulated system. As a result Tg and 
T^ , are much larger than the classical period, Tg , and Ts;^'^ . 
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The time modification factors Mjf'''^ and MJ^^^ are related as, 

mP = -3M^^^) = 3a, (81) 

wliere 

Hence, we conclude that the modification factors M^'^) and M'f^'> are directly 
proportional to the square of the nonlinearity, C^ in the system. From equa- 
tions (80) and (81) we note that as the quantum revival time T^ reduces 
by SaTo , T^ increases by aT^ . 

In the asymptotic limit, that is for (^ approaching zero, the quantum re- 
vival time in the modulated and unmodulated system are equal and infinite, 
that is T{ = Tg = oo. Furthermore, the classical period in the modu- 
lated and unmodulated cases are related as, T<^ = Tq A, as discussed in 
section IX.2.1. 

IX. 2. 3 Strong nonlinearity 

For relatively strong nonlinearity, the classical period, T^^ , and the quantum 
revival time, T^ , in the modulated system are related with the, Tq , and, 
Tq , of the unmodulated system as. 

The time modification factors, MJ^^^ and Mjf^\ follow the relation 

M^^) = M^^') = -/3, (84) 

70 



where 

Hence, for a relatively strong nonlinear case, j3 may approach to zero. 
Thus, the time modification factors vanish both in the classical and the 
quantum domain. As a result, the equations (75) and (76) reduce to, T^ = 
Tq A, and T^ = Tg , which proves the equality given in equation (83). 

The quantity (3, which determines the modification both in the classical 
period and in the quantum revival time, is inversely depending on fourth 
power of Planck's constant h. Hence, for highly quantum mechanical cases, 
we find that the recurrence times remain unchanged. 

X Dynamical recurrences in the Fermi accel- 
erator 

For the atomic dynamics in the Fermi accelerator, the classical frequency, u, 
and the nonlinearity, (, read as 

--(art) (««" 

and 

/ TT \ 2/3 

respectively. 

In the undriven system, we define the classical period as Tg = 2it/uj, 
and the quantum revival time as Tq = 2tx /JiC^. Hence, we find an increase 
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in the classical period and in the quantum revival time of the undriven Fermi 
accelerator as the mean quantum number r increases. 

The calculation for the matrix element, V, for large A^, yields the result, 

^'^-^. (88) 

and the resonance, N, takes the value 



N = ^^. (89) 

TT 

The time required to sweep one classical period and the time required for 
a quantum revival are modified in the driven system. In the classical domain 
and in the quantum domain the modification term becomes, 

and 

^ Aen) (1-/XT ^ ^ 

respectively, where 



Therefore, the present approach is valid for a smaller value of the modulation 
strength, A and a larger value of the resonance, N . 

A comparison between the quantum revival time, calculated with the help 
of equation (91), with those obtained numerically for an atomic wave-packet 
bouncing in a Fermi accelerator shows a good agreement, and thus supports 
theoretical predictions (Saif 2000e, 2002). 
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Figure 13: The change in revival phenomena for the wave packet originat- 
ing from the center of a resonance: (above) The square of auto-correlation 
function, C^ = | < ip{0)\ilj(t) > p, is plotted as a function of time, t, for an 
atomic wave packet. The wave packet initially propagates from the center 
of a resonance in atom optics Fermi Accelerator and its mean position and 
mean momentum are Zq = 14.5 and po = 1-45, respectively. Thick line corre- 
sponds to numerically obtained result for the Fermi accelerator and dashed 
line indicates quantum revivals for a wave packet in harmonic oscillator. The 
classical period is calculated to be An for A = 0.3. (bellow) The square of the 
auto-correlation function as a function of time for A = 0. As the modulation 
strength A vanishes the evolution of the material wave packet changes com- 
pletely. The wave packet experiences collapse after many classical periods 
and, later, displays quantum revivals at Tq . The inset displays the short 
time evolution of the wave packet comprising many classical periods for in 
the absence of external modulation. (Saif 2005a). 

The quantum recurrences change drastically at the onset of a periodic 
modulation. However the change depends upon the initial location of the 
wave packet in the phase space. In order to emphasis the fascinating fea- 
ture of quantum chaos, we note that for an initially propagated atomic wave 
packet in the nndriven Fermi accelerator the quantum revivals exist. How- 
ever, they disappear completely in the presence of modulation if the atom 
originates from the center of a primary resonance. In the present situation 
the atomic wave packet reappears after a classical period only (Saif 2005a), 
as shown in figure 13. 

We may understand this interesting property as we note that at a reso- 
nance the dynamics is effectively described by an oscillator Hamiltonian (Saif 
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2000c), 

Therefore, when (/? ^ 1, the Hamiltonian reduces to the Hamiltonion for a 



^ = -7rT + ^oCosv?. (93) 



harmonic oscillator, given as 

The Hamiltonian controls the evolution of an atomic wave packet placed 
at the center of a resonance. This analogy provides us an evidence that if 
an atomic wave packet is placed initially at the center of a resonance, it 
will always observe revivals after each classical period only, as in case of a 
harmonic oscillator. 

In addition, this analogy provides information about the level spacing and 
level dynamics at the center of a resonance in the Fermi accelerator as well. 
Since in case of harmonic oscillator the spacing between successive levels is 
always equal, we conclude that the spacing between quasi-energy levels is 
equal at the center of resonance in a periodically driven system. 

For a detailed study on the effect of initial conditions on the existence 
of the quantum recurrences, the information about the underlying Floquet 
quasi-energies, and the use of quantum recurrences to understand different 
dynamical regimes we refer to the references (Saif 2000c, 2005, 2005a). 
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XI Chaos in complex atom optics systems 

At the dawn of the 20*^* century, Max Planck provided the conceptual under- 
standing of the radiation fields as particles. The idea of the wave nature for 
material particles, coined by De Broglie, thus gave birth to the wave-particle 
duality principle. 

In later years the wave-particle duality became a prime focus of interest in 
the electron optics and neutron optics. However, easily controllable momenta 
and rich internal structure of atoms, provided inherent advantages to atom 
optics over electron optics and neutron optics. 

The enormous work of last three decades in atom optics has opened up 
newer avenues of research. These include many theoretical and experimental 
wonders, such as cooling of atoms to micro kelvin and nano kelvin scales, 
Bose Einstein condensation, and atom lasers. 

The study of the classical and the quantum chaos in atom optics started 
after a proposal by Graham, Shloutman, and Zoller. They suggested that an 
atom passing through a phase modulated standing light wave may exhibit 
chaos as the atom experiences a random transfer of momentum. 

Later, numerous periodically driven systems were studied in the frame 
work of atom optics to explore characteristics of the classical and the quan- 
tum chaos. On the one hand, study of these systems enabled the researchers 
to understand the theory of chaos in laboratory experiments, and on the 
other hand has brought to light newer generic properties, such as quantum 
dynamical recurrence phenomena in chaotic systems. 
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XI. 1 Atom in a periodically modulated standing wave 
field 

In the original Graham Schloutman and Zoller (GSZ) proposal, an altra cold 
atom in a phase modulated standing wave field exhibits Anderson like local- 
ization phenomenon in the momentum space. The exponential localization 
is visible as the momentum distribution exhibits exponential decay in the 
stochastic region. 

There exist domains on the modulation strength, where the GSZ system 
is almost integrable and the classical dynamics is regular. The exponential 
localization of the momentum distribution vanishes for these modulation 
strengths and occurs otherwise. Thus, the dynamical localization occurs over 
windows on modulation strength, separated by the domains which support 
integrable dynamics. 

The atom in the phase modulated standing wave has one more interesting 
feature. The classical phase space is a mixed phase space containing islands 
of stability in a stochastic sea. The modulation amplitude determines the 
ratio between the phase space areas corresponding to stability and chaos. 
This ratio serves as a measure of chaos. 

These dynamical features are experimentally verified by the group of 
Mark Raizen at Texas (Raizen 1999). In the experiment, sodium atoms 
were stored and laser-cooled in a magneto-optical trap. After turning off 
the trapping fields and switching on a standing light wave which has phase 
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modulation the atoms experience momentum kicks. The experiment demon- 
strates the exponential localization of the atoms in the momentum space 
(Latka 1995). 

The localization windows of GSZ model have a remarkable difference from 
the localization window of the Fermi accelerator, discussed in section VII. 2. 
In GSZ model, outside the localization windows, the dynamics is stable and 
bounded both classically and quantum mechanically. In contrast, in the 
Fermi accelerator the dynamics is bounded for modulation strengths below 
the localization window, whereas it is unbounded and chaotic for modulation 
strengths above the upper bound of the localization window both in the 
classical and quantum domain. 

XI. 2 Delta kicked rotor 

A rotor is a particle connected by a massless rod to a fixed point. While 
the particle rotates frictionless around this point it experiences an instanta- 
neous periodic force. The kicked rotor model has extensively been studied in 
classical chaos and at the same time it has become the beauty of quantum 
chaos. 

Historically the system provided the first study of dynamical localization 
in a periodically driven time dependent system (Casati 1979). Later, a formal 
connection was established between the kicked rotor model and one dimen- 
sion tight binding Anderson model with a time- dependent pseudo-random 
potential (Fishman 1982, Fishman 1984). This led to recognition that the 
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quantum suppression of classically chaotic diffusion of momentum of rotor is 
a dynamical version of the Anderson localization. 

The work of Mark Raizen, as discussed in section V.3, built a bridge be- 
tween atom optics and delta kicked rotor model (Raizen 1999). The work has 
led to the experimental verification of dynamical localization in delta kicked 
rotor (Moore 1995), to study the classical anomalous diffusion in quantum 
domain (Klappauf 1998), and to develop newer methods for cooling atoms 
(Ammann 1997). 

XI. 3 Ion in a Paul trap 

The idea of extension of two dimensional focusing of charged and neutral 
particles to three dimensions, led the invention of sophisticated traps and 
won the inventor, Wolfgang Paul, Nobel prize in 1989 (Paul 1990). A full 
quantum treatment of an ion in a Paul trap has been given in reference 
(Glauber 1992, Stenholm 1992). 

The trapped ions have been considered to generate the non classical states 
of motion of ions (Meekhof 1996, Monroe 1996, Cirac 1996), quantum logic 
gates (Monroe 1995, Cirac 1995), tomography (Leibfried 1996, Wallentowitz 
1995, Poyatos 1996, D'Helon 1996) and endoscopy (Bardroff 1996) of the 
density matrix . 

Two or more than two ions in a Paul trap may display chaos due to ion- 
ion Coulomb repulsion (Hoffnagle et al. 1988, Bliimel et al. 1989, Brewer 
et al. 1990). However a single ion moving in a trap in the presence of a 
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standing light field, displays chaos as well (Chacon and Cirac 1995, Ghaffar 
et al. 1997). The effective Hamiltonian of the system becomes 

where W{t) = a + 2qcost and Vq is the effective coupling constant. The 
parameters a and q denote the dc and ac voltage applied to the trap (Paul 
1990). 

The ion trap has been proposed to study the dynamical localization both 
in the position space and in the momentum space (Ghafar 1997, Kim 1999). 
In addition, the quantum revivals have also been predicted in the system 
(Breslin 1997). 

XI. 4 Fermi-Ulam Accelerator 

In the Fermi-Ulam accelerator model a particle bounces off a periodically 
vibrating wall and returns to it after observing reflection from a static wall, 
parallel to the vibrating wall. Based on the original idea of Fermi (Fermi 
1949), the classical Fermi-Ulam accelerator has been investigated extensively 
(Lichtenberg 1983, 1992). The system exhibits chaotic, transitional and pe- 
riodic solutions (Zaslavskii 1964, Brahic 1971, Lieberman 1972, Lichtenberg 
1980). 

The classical dynamics of the Fermi-Ulam accelerator model is expressible 
by a highly nonlinear mapping which connects the phase space position of 
the system at nth collision with the (n + l)th collision. The study of the 
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mapping revels that the orbit of a particle in the system can be stochastic 
(Zaslavskii 1964). However, the energy of the particle remains bounded and 
the unlimited acceleration as proposed by Enrico Fermi is absent (Zaslavskii 
1964, Pustilnikov 1983, 1987). 

In 1986, the quantum Fermi-Ulam accelerator was introduced (Jose 1986, 
Reichl 1986) which was later studied extensively to explore the quantum 
characteristics (Visscher 1987, Jose 1991, Scheininger 1991, Chu 1992). The 
beautiful numerical work of Jose and Cordery has demonstrated a clear tran- 
sition from Poisson statistics of quasienergies in the accelerator system to 
Wigner distribution (Jose 1986). The Poisson statistics implies integrability 
and Wigner distribution corresponds to nonintegrability (Brody 1981), hence, 
the transition is related to a transition from regular to irregular behavior. 

The study of quasi-energy spectrum of the quantum Fermi-Ulam system 
has led to establish that in general the time evolution of a quantum particle in 
the system is recurrent, and thus the energy growth remains bounded. How- 
ever, for some particular wall oscillations the particle may gain unbounded 
acceleration (Seba 1990). 

XII Dynamical effects and Decoherence 

Quantum interference effects are extremely fragile (Giulini 1996), and sen- 
sitive to dissipation and decoherence (Walls 1994, Zurek 1991). The phe- 
nomenon of dynamical localization is an interference effect. It is therefore 
extremely susceptible to noise and decoherence (Zurek 1995). 
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Figure 14: Comparison of the momentum distribution is made in the kicked 
rotor for the cases of (a) no noise, (b) 62.5% amphtude noise, and (c) dissi- 
pation from 13%/kick spontaneous scattering probabihty. Time steps shown 
are zero kicks (hght sohd hne), 17 kicks (dashed-dotted), 34 kicks (dashed), 
51 kicks (dotted), and 58 kicks heavy sohd for (a) and (c), and 0, 16, 32, 52, 
and 68 kicks, respectively, for (b). (Klappauf et al. 1998). 

First experiments (Bliimel 1996) using the dynamical localization of Ry- 
dberg atoms in microwave fields have clearly verified the destructive nature 
of noise. In atom optics the experimental realization of the kicked rotor has 
opened up a new door to study the influence of noise and dissipation on such 
quantum interference effects (Klappauf 1998). 

The presence of amplitude noise in kicked rotor destroys the exponential 
localization and broad shoulders develop in the momentum distributions. 
Furthermore, by tuning the laser field close to the resonance of the atom, 
the atom gets a non zero transition probability to jump to the excited state, 
and to spontaneously decay to the ground state. This spontaneous emission 
event destroys the coherence and thus the localization, as shown in figure 14. 

For more experiments on decoherence in dynamical localization we refer 
to (Ammann 1997). Moreover, there exists an extensive literature on the 
influence of spontaneous emission on dynamical localization, see for example 
(Graham 1996, Goetsch 1996, Dyrting 1996, Riedel 1999). 
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